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The problem of whether the electric charge is quantized or not has been longly 
investigated, both from experimental and theoretical points of view. Indeed, 
there are several experimental observations and theoretical arguments suggest- 
ing that the electric charges of the elementary particles we know are integer 
multiples of a fundamental unity: the charge of d-quark (— l/3e). 

Nevertheless, the existence of other particles carrying smaller electric charges, 
which electromagnetically interact with the ordinary photon, has not been ruled 
out so far; since these particles were supposed to have charge of order 10~^ e, 
thay are usually referred to as millicharged particles, or MCPs. 

In 1966 Kobzarev, Okun and Pomeranchuk proposed that there may be a 
mirror sector, that is a second particle sector identical to our one, but where 
weak interactions are right-handed; such a sector would provide to restore the 
strange lack of the left-right symmetry observed in nature. 20 years later, in 
1986, Holdom proposed a mechanism to induce MCPs via the kinetic mixing of 
two photons; this mechanism naturally applies to theories containing a mirror 
sector. 

If existing, mirror MCPs would have effects on the cosmological evolution of 
the universe (big bang nucleosynthesis, structure formation, cooling of stars and 
so on) as well on laboratory experiments, for instance the laser ones or those 
searching for dark matter. An other phenomenon where mirror MCPs may lead 
to observable effects is the positronium decay, as pointed out by Glashow in 
1986. 

The aim of these thesis is studying in detail the big bang nucleosynthesis 
(BBN) in models containing the mirror sector, which can be or not millicharged. 

The first part of this work (chapters [T] and [2]) consists in a review of basical 
concepts and formalisms used in particle physics and cosmology. Chapter [1] 
contains an introduction to the standard model of elementary particles, while 
chapter [2] is devoted to the standard model of cosmology. 

The problem concerning the quantization of the electric charge in theoreti- 
cal physics is dealt in chapter [3l three of the most challenging physics models 
(quantum standard model, grand unified theories and magnetic monopoles) are 
analyzed together with their predictions about the elementary particle electric 
charges; these charges end up to be quantized in all the three models. 

In chapter m the photon kinetic mixing, which leads the mirror matter to be 
millicharged, is introduced; then a review of the mirror theory and its cosmo- 
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logical features is given. We will argue that mirror MCPs are compatible with 
the theoretical bounds analyzed in the previous chapter. 

The last two chapters are devoted to the original work which has been done 
for this thesis. In chapterObig bang nucleosynthesis in the mirror scenario is an- 
alyzed in detail and the primordial abundances of light elements are worked out, 
but without considering at the moment the possible existence of millicharged 
interactions. 

Finally, in chapter El BBN in presence of the millicharged mirror sector is 
analyzed. After a review of the main processes leading to particle exchanges 
between the two sectors, we work out two approximated bounds on the photon 
mixing parameter e, which is related to the magnitude of the millicharges. Fi- 
nally, we propose a scenario worth to be analyzed in detail in future researches. 



Chapter 1 



The standard model of 
elementary particles 

The standard model of elementary particles (SM) is a field theory which de- 
scribes electromagnetic, weak and strong forces by mean of the group SU{3)c x 
SU {2)lxU and contains the fundamental particles which make up all mat- 
ter. 

In the literature the term "standard model" sometimes refers to the only 
electroweak sector, gauged by the group SIA{2)l x W(l)y; nevertheless, if not 
differently specified, by SM we will mean the whole SU{?>)c x SU{2)l x U{1)y- 
Its building required great efforts and contributes by various authors [1-3] , who 
received the Nobel Prize in Physics in 1979. 

The SM had several experimental successes; among them, it predicted the 
existence of several particles unknown at that time which were afterwards exper- 
imentally observed, like the 3 massive bosons which mediate weak interactions 
- and Z° - gluons, charm and top quarks. 

Nevertheless, the SM can not be considered a complete theory of fundamen- 
tal interactions because it leaves many unanswered questions. For example, it 
does not include gravity, which is the fourth known fundamental force, and it 
contains a wide numbe; of arbitrary parameter^ , which are chosen to fit the 
experimental data but can not be derived from first principles. 

In this section we briefly review the SM main features and its structure as 
a gauge theory. Wider treatments of these topics can be found in books, such 
as [4]. 

^There are at least 19: 3 gauge couplings - one possible choice is aem, ctg and sin^Ow 
-, 9 coupling constants of the Higgs boson with leptons and quarks, from which the fermion 
masses arise - see jjl.Sl 4 - three mixing angles and one phase - from the CKM matrix, 2 from 
the Higgs potential - one possible choice are the VEV v and the quartic coupling A - and 
finally the QCD 9 parameter. This number is still higher if neutrinos are massive particles, 
as indicated by several recent experimental and theoretical hints. 
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§J. Jl Symmetries and conservation laws 4 

1.1 Symmetries and conservation laws 

Let us consider a system specified by mean of several fields (f>j{x), j = 1, 
and their gradients d^(j)j{x); the theories we are considering can be derived 
applying a variational principle to the lagrangian density C, which is a function 
of and d^(l)j{x)\ 

£(x,t) = /:(0,(x),9^0,(x)) (1.1) 

The action integral S{VL), where VL is an arbitrary region of the four- dimensional 
space-time, is defined as: 

S{Q)= L{t)dt= f d^xC{(j)j{x),d^(j)j{x)) (1.2) 

The variational principle we will use is analogous to Hamilton's principle in 
mechanics and states that for any variations of the fields vanishing on the region 
surface, that is for any 

(pj{x) 4>j{x) + d(f)j{x), d(f)j{x G Q surface) = (1.3) 

the action has a stationary value: 

dS{Q) = (1.4) 

From this principle we can derive the Euler-Lagrange equations of motion: 

dC d dC _ 

In classic theories §| Noether's theorem holds [5]. It states that a conserved 
current defined as follows: 

J' = ^^^a (1.6) 

is associated with each generator of a continuous symmetry in the lagrangian. 
Moreover, the charge defined as: 

Q{t) = I £xJoix) (1.7) 

is a constant of the motion: 

Derivations of the Noether theorem can be found in any book about quantum 
field theory - see for instance [29-31]. 



^Generalization to quantum theories is not trivial, see §3.11 
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1.2 Gauge symmetries 

In the previous section we stressed the hnk of symmetries and conservation 
laws. We will show now that if the symmetry transformations are space-time 
dependent, they can be used to generate dynamics, that is, interactions between 
particles. These local symmetries are also called gauge symmetries. 

1.2.1 Gauge theories for abelian groups 

A group is called abelian if all the generators commute with each other. The 
simplest abelian group is the one- dimensional W(l), having the unity matrix 
as generator. Because of the important physical applications of this group - it 
enters both QED and SM - we will refer to the only U{1) in this section. 

If a system is invariant under the fields can be changed by an arbitrary 
phase 9: 

^ (j)'j{x) = exp-'*^ (j)j{x) (1.9) 

without any effects on the physically measurable quantities. In the equation 
above qj is the charge of the particle represented by the field in units of the 
electron's charge e. 

If 6 does not depend on the position in the four- dimensional space-time, the 
theory is globally symmetric; if 6 depends on x instead, the theory is locally 
invariant and is called gauge theory. In this case the transformation 

(f)j{x) (f)'j{x) = exp-'«^^(^) (f)j{x) (1.10) 

leaves C unchanged. 

Let us consider a lagrangian density C{(f),d^(j)) invariant under the global 
tranformation (11.91) . The terms in C containing only the fields are also invari- 
ant under the local transformation (11.101) . but the ones involving gradients are 
not because 

d^(f)j{x) exp~**^(^) d^(j)j{x) + ...^ exp-**^(^) d^(f)j{x) (1.11) 

To generalize a global invariance to a local one therefore, the derivatives must 
be replaced by special combinations called covariant derivatives which are 
chosen such that 

D^<j)^{x) exp-**^(^) D^<f)^{x) (1.12) 

The derivative L*^ depends on a field which is called gauge boson and is used to 
reabsorbe the extra terms in (II. lip . 
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QED 

In QED, which is gauged by the group W(l), the covariant derivative has the 
form: 

D^^ = d^ + teq,A^ (1.13) 
and the field A^, which is the vector potential, transforms as 

A^^A'^ = A^ + ^d^e (1.14) 

If we want the theory to be gauge invariant, we have to replace the simple 
derivative with the covariant one in the Dirac lagrangian for free spinors 

iplip— m)ilj = — eqjJjL — m)il) (1-15) 

This way we get interactions between fermions and photons through the extra 
term —eqjipj^ip. This phenomenon is also present in theories gauged by more 
complicated groups - see §1.2.21 in general transforming a global symmetry to 
a local one in a lagrangian for free fields , we get interactions between the fields 
and the gauge bosons. 

1.2.2 Gauge theories for non-abelian groups 

Let us now consider a non-abelian gauge group G having a set of generators Tj 
with j = 1, No; the TjS obey the commutation relations: 

[Ti,Tj] = ifijkTk (1.16) 

where fijk are called the structure constants of the group and are antisymmetric 
under interchange of any pair of indices. 

If the fields tranform according to some representations of G, the generators 
Tj will be represented hj nxn matrices Lj; the field transformations are specified 
by Nc parameters which we will call 6j and can be written as: 

0^ 0' = exp-^^'V (1-17) 
where is a multiplet of fields 



/ 01 \ 

02 



1.18) 
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Similarly to the abelian case, the lagrangian can be invariant under a global 
transformation U{6j). But, if the parameters 9j depend on the position x, the 
gradients contained in the lagrangian will break the invariance at local level. 
This problem is avoided defining a covariant derivative with the property 

D,cl>{x) ^ D'^ct>'{x) = U{9)D^ct>{x) (1.19) 

Assuming that the lagrangian contains gradients only through the covariant 
derivative, the invariance under local non-abelian gauge transformations is en- 
sured. Generalizing the procedure we adopted for QED, we can introduce one 
vector field VF^ (x) for each generator of the group and then write the covariant 
derivative as 

= [d^ + igl.- W^(x)]0(x) (1.20) 

where g plays the role of a coupling constant. Finally, if we impose that the 
property fll.lQp must be vahd also for defined in (11.201) . we will derive how 
the gauge fields W^ (a;) must transform: 



L ■ w = u{e) 



L-W, + -U-\9)d^i9)Ui9) 
9 



U-\9) (1.21) 



As disclosed at the end of §1.2. H passing from a global to a local symmetry 
induces interactions in free particle lagrangians for both abelian and non-abelian 
gauge theories. The SM is based on this formalism, as we will see in the following 
sections. 



1.3 Spontaneous symmetry breaking 

Spontaneous symmetry breaking (SSB) occours when a lagrangian has some 
exact symmetry but the solutions of the problem do not. The Goldstone 
theorem states that the spontaneous breakdown of a continuous symmetry in 
field theory implies the existence of massless spinless particles. These particles 
are referred to as Nambu- Goldstone bosons or simply Goldstone bosons. Their 
number is equal to the number of the broken generators. The Goldstone theorem 
has been first studied by Nambu [6-8], and later proved by Goldstone and 
others [9,10]. 

Dynamical systems in which the ground state does not possess the same 
symmetry properties of the lagrangian are very interesting for particle physics. 
Indeed it comes out that when SSB occours in a gauge theory involving massless 
vector fields and scalar fields, the Goldstone bosons disappear and re-emerge as 
the longitudinal mode of the vector fields, which therefore behave like massive 
particles {Higgs phenomenon) . This way we get the three massive bosons - 
are - which mediate weak interactions. 
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1.3.1 SSB of a global symmetry and Goldstone theorem 

In this section an example of SSB in classical field theory is briefiy reported. 
Let us consider a complex scalar field with the lagrangian: 

CssB = {dM9^<l^*) - - (1-22) 

which is invariant under the global transformation 

(f)'{x) = exp-'^ (1.23) 

£ can be thought as the sum of a kinetic term and a "potential" one: 

V{(f)) = fi^4>^* + A(00*)2 = fi^p + Ap^ p = (1.24) 

Since the kinetic term in L vanishes for = const, minimum points for V are 
also minima of the total energy and are therefore called vacua or ground states. 
Clearly V has a minimum only if A > 0. Assuming A > and p^ > 0, V has one 
minimum in p = 0. If /i^ < instead, that is p does not represent a physical 
mass for the field 0, there are infinitely many minimum points for V in 



(j)vac = -^exp'^ where ^ ^ y ~X~ (1-25) 

and A is a real arbitrary number. All the minima are equivalent with each other, 
therefore anyone of them can be arbitratly chosen. In particular we can choose 
the one laying on the real axis, having A = 0, and therefore get: 

(f){x) = ^[v + ax)+txix)] (1.26) 

Sobstituting the field in ( 11.26^ in the lagrangian in ( 11.22^ we get, ignoring some 
constant terms: 



1 



>C = 7;id,0' + 7;id,xr - At^T - Xv^e + x') - -K^ + xY (1-27) 



This way we rearranged £ as a functions of two fields x ^md ^, where ^ is 
massive: 

mj = 2Xv^ (1.28) 
The massless field x is the Goldstone boson. 
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1.3.2 SSB of a local symmetry and Higgs phenomenon 

We saw in §1.2.11 and 11.2.21 that the imposition of a local symmetry implies the 
existence of a certain number of massless gauge bosons. If we want these bosons 
to be massive and introduce this feature explicitly by a symmetry-breaking mass 
term, we will spoil the renormalizability of the theory. The high momentum 
limit of the propagators indeed is dominated by 



which leads to divergent terms. 

A solution to this problem is found spontaneously breaking the involved 
symmetry. This way we get a certain number of Goldstone bosons, which we 
would like to eliminate. Remarkably the Goldstone theorem is evaded in gauge 
theories because its derivations requires axioms which are not all compatible 
with the gauge-fixing condition. 

What finally comes out is the disappearance of the Goldstone bosons, which 
re-emerge as longitudinal DOFs of the gauge bosons, which therefore becomes 
massive. 

This phenomenon was first studied by Anderson [11] for condensed- matter 
physics and later generalized to relativistic field theory by Englert and Brout [12] 
and by Guralnik et al. [13]. But it is known as Higgs phenomenon because this 
author gave the most complete treatment [14,15]. The first demonstration that 
field theories are renormalizable also in the presence of SSB has been provided 
by 't Hooft [16]. 

A short example 

Let us consider again the lagrangian density fll.22p : local invariance is achieved 
substituting ordinary derivatives with covariant ones; moreover, a kinetic term 
for the gauge bosons must be added. This way, the W(l) gauge invariant la- 
grangian becomes: 



This lagrangian has in total four DOFs, two from the massless vector boson 
and two from the scalar complex field. As for the globally invariant case, if we 
assume A > and fi^ < there is a ring of degenerate ground states. If we 
re-defyne the scalar field as: 



^2 _ ^2 



const 



(1.29) 




- A(#*)^ 



(1.30) 





(1.31) 
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and substitute it in the lagrangian fll.30p we get: 

^ = -If.^F^" + + ^(^.0' + lid.X? - Xv'e + ... (1-32) 

which seems to describe the interactions of a massive vector field - having three 
DOFs - and two scalars - having a total of two DOFs. Therefore we have in 
total five DOFs and not four as in the original lagrangian in eq. fll.30p . This 
one more degree of freedom is only apparent and can be eliminated choosing a 
particular gauge called the U gauge in which the unphysical DOF is absorbed 
in the arbitrary field phase. 

Let us see how this mechanism works in a W(l) gauge theory. Since the 
parameter 9{x) can be arbitrarly chosen, we can set it equal to the phase of 
4>{x) at each space-time point, so that: 

= exp-*^(^) = ^[v + r]{x)] (1.33) 
V2 

This way both (p' and rj are real and the lagrangian becomes 

^ = -\f'^.F"^^ + '^A'^A'^ + \{d,vr - XvW + ■■■ (1.34) 

where 

a;(x) = A,{x) + -d,e{x) 

F',. = d,Al - d,A'^ (1.35) 

This way C describes the interactions of a massive vector boson A'^ and a real 
scalar field t], called the Higgs boson with mass 

ml = 2Xv^ = -2i? (1.36) 

In conclusion, when a symmetry is spontaneously broken the gauge boson ac- 
quires a mass while the Goldstone boson disappears, leaving the Higgs boson as 
the only present scalar field. The non-abelian extensions of this mechanism are 
called Yang-Mills theories and have a great importance in the SM construction. 

1.4 Parity violation in weak interactions 

Let us introduce the chirality projection operators: Pr/l'- 

Pr = ^(1 + 7') =^ Pr^ = ^r 
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Pl = ^(1-7') =^ Pii^^i^L (1.37) 

where ipR/L are called respectively right- and left-handed components of a spinor. 
It can be shown that for a massive fermion handedness does not commute with 
the Dirac hamiltonian: 

i^Dirac = (« ■ P + /?m) (1.38) 





and therefore it is not a good quantum number. 

In 1956 Lee and Yang proposed that weak interactions may be not invariant 
under parity transformations [17]. Actually, the problem of wether parity is 
conserved or not had been investigated in several contexts at that time but Lee 
and Yang where the first who observed that there were no experimental tests of 
parity conservation in weak interactions. Moreover, there was evidence for the 
decay in two modes {r — 9 paradox) in which the final states have opposite 
parities. 

In their paper they also proposed several possible experimental tests; the 
first one was carried out by Wu et al. [19], who studied P decay of ^°Co. They 
found a clear violation of parity conservation: weak interactions are left-handed. 

Let us now introduce the quantum number helicity, which is defined as: 



S • P = S ..p (1.40) 




Helicity is therefore a measure of the allignement of the spin and momentum 
vectors; the corresponding helicity projection operators are: 

n±(p) = 1 (1 ± s . p) (1.41) 

Helicity commutes with the Dirac hamiltonian (11.381) and is therefore a good 
quantum number for spinors. It can be shown that in the high energy limit, 
that is for m = 0, the helicity projection operators become: 

hmH±(p) = i(l±75) =P«/i (1.42) 

that is, they coincide with the chirality projection operators. 
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1.5 Construction of the SM 

This section's aim is to introduce some methods and ideas which will be useful 
later on in this thesis. It does not give an exhaustive nor complete introduction 
to the standard model, which can be found in textbooks. 

The electro-weak interactions are gauged by the group 5W(2) xW(l), having a 
total of 4 generators: Y, that is the unit matrix 1 divided by 2, and r^, i = 1,2, 3, 
that is the Pauli matrices cTj also divided by 2. 

The color group SU{3), when gauged, gives rise to quantum chromodynamics 
- or QCD, which describes the strong interactions of colored quarks and gluons. 
This group has 8 generators Tq, a = 1,...,8 which corresponds to 8 gauge 
fields Gi, i = 1, . . . , 8 which are called gluons. Since the color symmetry is 
unbroken, gluons do not acquire any masses. They do not take part in electro- 
weak interactions, that is they are singulets of the symmetry SU{2) x W(l), but 
they carry color. They can therefore interact with each other as well as with 
quarks, which carry color too, but not with leptons, which are color singulets. 

QCD has two peculiar properties; the first one, called Asymptotic freedom, 
was discovered in the early 1970s by David Politzer [20] and by Frank Wilczek 
and David Gross [21]. Asymptotic freedom means that in very high-energy 
reactions, quarks and gluons interact very weakly; this implies that in high- 
energy scattering the quarks move within nucleons, such as the neutron and 
proton, essentially as free, non-interacting particles. 

The second property, called Confinement, implies that the force between col- 
ored particles (such as quarks) does not diminish as they are separated. Because 
of this, it would take an infinite amount of energy to separate two of these parti- 
cles, which thus can not be isolated and are confined in bound states of neutral 
net color. Although analytically unproven, confinement is widely believed to be 
true because it explains the consistent failure of free quark searches, and it is 
easy to demonstrate in lattice QCD. 

Let us consider now the only electro-weak part of the theory; in a model 
with four generators we have four gauge bosons, which we call: 

forW(l) 

Wl i = l,2,3 for5W(2) (1.43) 

If we want three of them to be massive, we need at least 4 independent scalar 
fields. The simplest choice of them {minimal standard model) is a douplet of 
complex scalar fields: 

^ = ( / ) (1-^*) 
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We want the model to include also leptons, which are electron, muon and tauon 
with the corresponding neutrinos and anti-particles, and quarks; neutrinos are 
massless §, that is they are helicity eigenstates. Moreover they experience only 
weak interactions, which are - see §1.41 - left-handed; therefore the only left- 
handed part of the neutrino has to enter the model and leptons are represented 
as 3 

^^=[''e) e« = (e)«~(l,l,-2) (1.45) 

The first is a SU{2) douplet containing left-handed neutrino and electron which 
are at the present both massless. The second is a singulet representing the right- 
handed electron, also massless. The brackets indicate the SU{3) and SU{2) 
content of the multiplets and their U{1) hypercharges. 
Similarly u and d quarks enter the model as: 



tiK = (ti)H ~ (3, 1, da=(d)a~ [3,1,-1^ (1.47) 

Moreover, leptons have a global lepton charge L = 1 and a family lepton charge 
Le = 1, while quarks have a baryon charge i? = |, so that baryons consisting 
of three quarks have B = 1. 

The charge conjugation matrix C makes the transformation 

fR = ci,fi rL = c^orR (1.48) 

where / stands for a generic fermion field - lepton or quark. The fields marked 
with have opposite gauge charges as well as opposite chiralities with respect 
to fermions. Also barion and lepton numbers are assigned with opposite signs. 
The tranformations of any fields / under the gauge groups are: 

iY^6(.) / lYw0{x)\ / 
f ^ f = exp 2 J ~ 1 _ 1 —-^ f (under W(l 



■^There is now convincing evidence that neutrinos change from one flavor to another, what 
imphes that neutrinos have non-zero masses. There are several possible ways to add neutrino 
masses to the model, for instance via a Dirac or a Majorana term. Nevertheless, theories 
containing such masses are considered beyond the standard model, which contains instead no 
mass terms for neutrinos. 

^There are three lepton families in the SM; nevertheless, since they have the same proper- 
ties with respect to the SU{3) x SU{2) x U{1) symmetries, we consider for sake of semplicity 
the only electronic family with the corresponding light quarks u and d. 
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f^f = exp-^^-^/~ (^1-z/H.^^^j/ (under 5W(2)) (1.49) 

where Yw and Iw are the field eigenvalues of the hypercharge Y and the weak 
isospin /. These numbers, together with the transformation properties under 
SU{2)l and SU{3)c and the electric charge are reported in Table ll.li We 



Table 1.1: The SM particle content. 
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analyze below the main interactions occourring in the model; most of them 
naturally arise when we impose gauge invariance - and therefore replace ordinary 
derivatives with covariant ones - in the lagrangian: 

C = CDiracik, Cr) + CMaxweu{B^, W^) + CsSb{(1)) (1-50) 

Coupling of the gauge bosons to the Higgs scalars: the electric charge 

Coupling of the gauge bosons (GB) to the Higgs scalars (S) are found imposing 
gauge invariance to the SSB lagrangian, which therefore becomes: 

(/.j^x (1.51) 

The potential involving the Higgs fields is the same than in §1.3.11 and its pa- 
rameters A and /z^ are chosen so that there are infinitely many minima: 

\/(0t0) = ^2^t0 + A(0V)^ A > 0, /i^ < (1.52) 

The most common choice for the ground state is: 

0™c = ^ ^ ^ (1-53) 



GB~S 



+ i -W„ -T + i -15„ 
M ' 2 ^ 2 
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None of the infinitesimal transformations in (11.491) leaves it unchanged and thus 
it seems that all the guage bosons acquire a mass via the Higgs mechanism. But 
we know that the photon is massless! Paying more attention we can see that 
the combination of generators 

Q = rsL + ^ (1.54) 

has no effects on (p^ac- To identify the massless vector boson - that is the photon 
- it is therefore useful to re-define the fields as: 

= cos 6wA^ + sin OwZfj, 

Wjl = smOwA^- cosOwZf, (1.55) 

where 9]y is the Weinberg angle, which is a free parameter with value sin^ 9\y ~ 
0.223 chosen to fit the experimental data. Making the calculations we can see 
that couples to the only Q if 

e = g sin 9w = g' cos 9w (1.56) 

With these substitutions in the lagrangian in (11.511) we finally get a theory in 
which: 

• A^ remains massless because it couples to only through the unbroken 
generator Q. 

• Q measures the electric charge in units of e. 

• The two charged generators defined as 

= {W^ T W^) (1.57) 
and the neutral one becomes massive via the Higgs mechanism with 

Mw = ^-^ Mz = ,r^ (1.58) 

2 2 cos 9w 

• The Higgs scalar which survives becomes massive with m^/ = \/ —2fj,'^. 
Coupling of the gauge bosons to the leptons 

Interactions of the gauge bosons (GB) and the leptons (1) arise when imposing 
gauge invariance to the Dirac lagrangian for massless spinors: 

Cgb-i = IlII^ [d^ - ^ig'B^ + hgr ■ h + em^^^d^ - ig'B^)eR (1.59) 
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Left-handed components of the leptons couple to charged bosons W . 



;i.60) 



while both left- and right-handed components couple to the neutral bosons 

Ci-A = -g sin 9w {eil^ei + eR-f'^eR) A,, = -eJ'^^j^fA^ (1.61) 



and Zl: 



i-z 



e tan 6' VI/ 

e 



csc^ 9w {^Ll^^L - glI^cl) + e7^e 



^3 - 2sin^ OwJqed 



2 cos 9w sin 9w 
where we used the electro-magnetic current 

and the third component of a triplet of weak-isospin currents: 

^3^ = hrcr^h 



;i.62) 



1.63) 



1.64) 



Coupling of the fermions to the Higgs: fermion masses and CP vio- 
lation 

The last interaction we analyze arises from the Yukawa term, which must be 
added by hand. Let us now take in account of all the three fermion generations; 
the Yukawa term has the form: 



Yukawa 



-Y:^ {hAeR^j) - Y^' {qL,k<l>dji,i) - Y^ {qL,n.<pUR,n) + h.c. 



kj = d,s,b 



m,n = u,cA 



1.65) 



Diagonalizing the three matrices in (11.651) we get 9 of the SM free parameters: 

(1.66) 



D 

e,d,u 



( ye,d,u 

y^^,s,c 
V yr,b,t 



which are proportional to the fermion masses: assuming indeed = (p^ac we get 
for instance from the leptonic term in (11.651) a mass for the electron {fi, r): 



£vac 



Y,v_ 



Yukawa— el 



72 



ee 



m 



_ Y^ 



1.67) 
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without explicitely break the symmetry in the lagrangian by mean of a Dirac 
mass term. The same procedure apphes to quarks; note that, in this model, the 
neutrino remains massless. 

But we can not perform the diagonahzation of the matrices Y for free; indeed 
we must use four ausihary matrices (per generation) Ul/r and Vl//? such that: 

y,^ = UlY.Un = VlY.Vn (1.68) 

but the three quark fields (per generation) we have are not enough to re-absorbe 
all them; one of them survives and is responsible for the flavor violating charged 
currents [Cabibbo-Kobayashi-Maskawa or CKM matrix). 

This matrix leads the SM to have 4 more physical parameters, one of which 
is a phase which makes the CKM matrix not real. It can be shown that the CP 
invariance, which is valid in the leptonic sector of the SM, does not hold in the 
hadronic sector because the presence of this phase of the CKM matrix. A small 
CP violating effect was first observed in 1964 studying the and decay - 
see [22] for a more detailed review. 



Chapter 2 

Standard cosmology 



In this chapter some topics of standard cosmology which will be used in chapters 
[5] and [6] are summarised. These topics can be also found in many textbooks and 
reviews - see for instance [75,76,83]. 

2.1 The expanding universe 

The Friedmann-Robertson- Walker metric 

The cosmological principle states that on large scales the Universe is to a good 
approximation homogeneous and isotropic. By homogeneity we mean that the 
geometrical properties are the same at all spatial locations, while by isotropy 
we mean that the geometrical properties do not sigle out any special direction 
in space. 

Assumed that the cosmological principle holds, it is convenient to describe 
the universe using a particular coordinate system, appropried to the special 
class of observers for whom the universe appears isotropic. 

The most general space-time metric describing a Universe compatible with 
the cosmological principle can be determined entirely by symmetry considera- 
tions without any references to the energy sources or Einstein's equations. It is 
called Friedmann-Robertson- Walker metric (FRW) and can be written as: 



where {t, r, 9, (p) are comoving coordinates □, a{t) is the cosmic scale factor, 
and, with an appropriate rescaling of the coordinates, k can be chosen to be 
—1, or for spaces of constant positive, negative, or zero spatial curvature, 

^ A "comoving observer' follows the expansion of the Universe including the effects of any 
inhomogeneities that may be present. 



ds^ = dt^ - a^{t) 



dr"^ 





1 — fcr^ 
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respectively. The coordinate r is dimensionless, while a{t) has dimensions of 
length. The conventions for signs we used are reported in appendix I A. 1[ 

The value of k and the form of the function a{t) can not be determined by 
geometrical consideration; they have to be determined with Einstein's equations 
once the matter distribution is specified. 



Redshift 

Let us consider an expanding universe. The light emitted by a distant object 
is described in the quantum mechanics scenario in terms of freely-propagating 
photons. Since the wavelength of a photon is 

A = - (2.2) 
P 

and the momentum p changes in proportion to a^^, the wavelength at time t^, 
denoted as Aq, will differ from that at time t, denoted as A, by 

^ = 4^ (2-3) 
Ao a[to) 

As the Universe expands, the wavelength of a freely-propagating photon in- 
creases, just as all physical distances increase with the expansion. 

Hence, we introduce a new variable related to the scale factor a which is 
directly observable, the redshift of an object, z: 

z = ^^ l + z = ^ (2.4) 

A a[t) 

where Aq is the detected wavelength and A is emitted wavelength. Any increase 
(decrease) in a(t) leads to a red shift (blue shift) of the light from distant sources. 
Since today astronomers observe distant galaxies to have red shifted spectra, 
we can conclude that the Universe is expanding. 



The Hubble's law 

The Hubble's law is a linear relationship between the distance of an object and 
its observed red shift; it can be expressed as 

z^ Hd^ lO^^^cm^i X d (2.5) 

where d is the proper distance^ of a source, and H is the Hubble constant or, 
more accurately, the Hubble parameter (because it is not constant in time, and 

^The proper distance of a point P from another point Pq is the distance measured by a 
chain of rulers held by observers which connect P to Pq- 
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in general varies as t defined by 

d(t) 

if = 4t 2.6 

a{t) ^ ' 

At present time the Hubble parameter value Hq is not known with great accu- 
racy, so it is indicated by 

Hq = 100 h km s-^Mpc"^ h = 0.73 ± 0.02 (2.7) 

where the dimensionless parameter h contains the uncertainty on Hq [38]. 

The present age and the local spatial scale for the Universe are set by the 
Hubble timen and radius 



Hq^ ~ 9.778 X lO^h-^ yr ~ 3000 Mpc (2.8) 

The Friedmann equations and the equation of state 

The expansion of the Universe is determined by the Einstein equations 

Gal3 = Rap — -RgaP = ~Taf3 + J^QaP (2.9) 

^ c 

where Rap is the Ricci tensor, R is the Ricci scalar, gap is the metric tensor. 
Tap is the energy- momentum tensor, and A is the cosmological constant 0. For 
the FRW metric (12. ip , they are reduced to the form (see appendix lA.ip 

h Atv 

- = --G(p + 3p) (2.10) 
a 3 

for the time-time component, and 



a 



+ 2(-) +2- = AtxG{p-p) (2.11) 
\aj a"^ 



for the space-space components, where, if the cosmological constant is present, 
p and p are modified according to flA.9p . From equations fl2.10l) and 02.111) we 
obtain also 

a J 3 



■^Note that earlier than some time, say tx, or better for a less than some ax, our knowledge 
of the Universe is uncertain, so that the time elapsed from a = to a = ax cannot be reliably 
calculated. However, this contribution to the age of the Universe is very small, and most of 
the time elapsed since a — Q accumulated during the most recent few Hubble times. 

^The cosmological constant A was introduced by Einstein for the need to have a static 
Universe. However, now we know for sure that the Universe is expanding. In fact, one of the 
biggest theoretical problems of the modern physics is to explain why the cosmological term 
is small and not order Mp. 
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Equations fl2.10p and fl2.12p are called the Priedmann equations; they are 
not independent since the second can be recovered from the first if the adiabatic 
expansion of the Universe is taken into account. The Friedmann equation (12.121) 
can be recast as 

^ ^ -l = n-l (2.13) 



where Q is the ratio of the density to the critical density necessary for closing 
the Universe § 

, Pc^l^ (2.14) 

Since H^a^ > 0, there is a correspondence between the sign of k, and the sign 
of {Q - 1) 

CLOSED k = +l =^ n>l 

FLAT k = =^ n = l 

OPEN k = -l =^ n<l 

From equation fl2.13p we find for the scale factor today 

which can be interpreted as the current radius of curvature of the Universe. If 
r^o = 1, then oo has no physical meaning and can be chosen arbitrarily (it will 
always cancel out when some physical quantity is computed). 

In order to derive the dynamical evolution of the scale factor a{t), it is 
necessary to specify the equation of state for the fluid, p = p{p). It is standard 
to assume the form 

p = wp (2.16) 

and consider different types of components by choosing different values for w. 

From equation (12.101) . models of the Universe made from fluids with — 1/3 < 
w < 1 have d always negative; then, because today a > 0, they possess a point 
in time where a vanishes and the density diverges; this instant is called the Big 
Bang singularity. Note that the expansion of the Universe described in the Big 



^The present value of the critical density is poc = 1-88 x 10^^^ g cm^'^; taking into 
account the range of permitted values for h, we get poc ~ (3 — 12) x 10~^^ kg m""^, which 
corresponds to a few H atoms/m'^. Just to compare, a 'really good' vacuum (in the laboratory) 
of 10~^ N/m^ at 300 K contains ^ 2 x 10^^ molecules/m'^. The Universe seems to be empty 
indeed! 
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Bang model is not due in any way to the effect of pressure, which always acts 
to decelerate the expansion, but is a result of the initial conditions describing a 
homogeneous and isotropic Universe. 

If the Universe is filled with pressureless non-relativistic matter (dust), we 
have p <ti p and thus w = 0. Instead, for radiation, the relativistic ideal gas 
equation of state p = l/3p will be used, from which follows w = 1/3. Another 
interesting equation of state is p = —p, corresponding to w = —1. This is the 
case of vacuum energy , which will be the relevant form of energy during the 
so-called inflationary epoch. 

The a = component of the conservation equation for the energy-momentum 
tensor, T°^^ = 0, gives the 1st law of thermodynamics: 

d{pa^) = -pd{a^) (2.17) 



where the change of energy in a comoving volume element is equal to minus the 
pressure times the change in the volume. 

From equations (12.161) and (I2.17p . assuming w independent of time, we can 
obtain the relation 



poca-=^(i+™) (2.18) 



Some special cases are: 



Radiation : Pr = ^Pr =^ Pr oc a 

Matter : pm = =^ pM oc 

Vacuum energy : Pv = —pv ^ oc const (2-19) 

The overall density in the Friedmann equation is: 

P = Pr + Pm + Pv (2.20) 

At the very beginning, the universe was radiation dominated; then it became 
matter dominated and, in absence of vacuum energy, it will continue to be 
matter dominated. 

2.2 Flat models 



From equations (12.121) and fl2.18p we get for /c = 



ao 



(2.21) 
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Now we shall find the solution to this equation appropriate to a Rat Universe. 
For fig = 1 integrating equation (12.211) gives 

a{t)=ao[^-J (2.22) 

which shows that the expansion of a flat Universe is indefinetely long in time; 
equation (I2.22p is equivalent to the relation between the cosmic time t and the 
redshift 

t = to(l + ^)"'^'^"'^^' (2.23) 
From equations (12.221) . (I2.23P and (12.181) . we can derive 

H ^ ^ = --^ = Ho'-^ = Hoil + zf'^-^/' , (2.24) 
a 3{1 + w)t t 

ad 1 + ?)W ,^ 
q = -— = — - — = const. = go , (2.25) 
2 

^0 = on ^ m ' (2.26) 

3(1 + w)Ho 



In appendix lA. 21 we report the above relations for the special cases of a Universe 
dominated by matter or radiation. A general property of fiat Universe models is 
that the scale factor a grows indefinitely with time, with constant deceleration 
parameter go- The pressure role can be stressed again by observing that increas- 
ing w, and therefore the pressure, leads also to an increase of the deceleration 
parameter. 

A cosmological model in which the Universe is empty of matter and has a 
positive cosmological constant is called the de Sitter Universe. From equations 
f [0|) and fimi) we obtain 

a\2 SttG 

= ^PA (2.28) 



which for positive A implies the exponentially fast expansion 

a(t) = aoe^(*-*o), ^ = ^/^a = f (^.29) 

corresponding to a Hubble parameter constant in time. In the de Sitter vacuum 
Universe test particles move away from each other because of the repulsive 
gravitational effect of the positive cosmological constant. 
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Finally, the age of a flat Universe containing both matter and positive vac- 



uum energy {Q = + = 1) is 



1 



In 



1 + 



,1/2 



PA 



(2.30) 



It is interesting to note that, unlike previous models, a Universe with > 0.74 
is older than Hq^ because the expansion rate accelerates. Also, when il\ 1 
the time to ^ oo. For this reason the problem of reconciling a young expansion 
age with other independent age determinations (like for example the globular 
clusters) has often led cosmologists to invoke a cosmological constant. 



2.3 Equilibrium thermodynamics 

In this section we will study the properties of the Universe considered as a ther- 
modynamic system composed by different species (electrons, photons, neutri- 
nos, nucleons, etc.) which, in the early phases, were to a good approximation in 
thermodynamic equilibrium, established through rapid interactions. Of course, 
going back to the past, the cosmic scale factor decreases while the temperature 
becomes higher. 

To work out the physical processes at some time t, we need the distribution 
function /^(x,p,t) of the present particle species. We assume, coherently with 
the hypothesis a homogeneous universe, that /a does not actually depend on 
the coordinates x, therefore /A(x,p,t) = /A(p,t). 

All interactions that involve elementary particles have a short range 0. We 
therefore assume that the role of these interactions is providing a mechanism 
for thermalization without affecting the form of the distribution function (ideal 
gas approximation) 0. We can therefore use the equilibrium Bose-Einstein or 
Fermi-Dirac distribution functions: 

/a(p, t)£p = ff^ ■ (2-31) 
l^^J exp^I^ ±1 

where gspinA is the spin-degeneracy factor of the species A, is the chemical 
potential, the signs + and - correspond respectively to fermions and bosons and 
Ep is given by the mass-shell equation: 

El = p2 + (2.32) 

^The only exception is the Coulomb force, that is anyway shielded in plasma by the Debye 
effect. 

^The ideal gas approximation is valid at low densities, that is p < where r is the 
average distance between two particles. 
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From the distribution function in (12.311) it is straightforward to calculate the 
number density n, the energy density p and the pressure p for every particle 
species in thermal equilibrium: 



J ZTT JniA exp — T — ihl 



irriA exp^^ ±1 



PA = I Ef{k)d'k = ^ I dE (2.34) 



r {E 



2 ^2 



ml 2 



= - y kt;(k)/(k)rf^k ='-^i dE (2.35) 



exp T ±1 



2.3.1 Entropy and energy 



Entropy per comoving volume can be defined, up to an additive constant, as 
(see ^KM- 

This quantity is conserved, that is 

dS = (2.37) 
We will assume below that p <^T and use the entropy density s, defined as: 

s{T) ^ ^ = %<lUT)T' (2.38) 

bosons fermions 

This way, the functions qf,/f{T), called entropy DOFs, contain all the integrals 
in (12.361) . This formalism is convenient because, as we can see from equations 
(IA.23p . the energy density and the pressure of a non relativistic species is expo- 
nentially smaller than that of a relativistic species and thus we can often ignore 
its contribution. 

An analogous formalism can be used for the total energy density p, which 
can be defined as: 

piT) = ^9tot{T)T' (2.40) 
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gUT) ^ E9,{T)(^)\l E 9fiT)(^V (2.41) 

bosons fermions 

The functions gh/f(T) are analogous to the functions qb/f(T) and are called 
energy DOFs. For ultrarelativistic species q(T) and g(T) are constants of the 
temperature and we have: 

QspiriA ~ 9A = QA 

g^ = 2.3 = 6 gj = 2 g^^ =2-2 = 4: (2.42) 

When photons, neutrinos and electron-positrons are in thermal equilibrium, 
that is T ^ 2 ^ 3MeV, the number of DOFs is easy to work out: 

gtot = qtot = 2 + 1- (4 + Q) = 10.75 (2.43) 

o 

Also, when the annihilation process is over, we get: 

Tj,/4\i 7 /4\t 

r=(n) = '^r'iu) 

= 2 + 5 ■ ■ (n) ~ ^-^^ 

The temperature factor (^-^^ ^ for neutrinos will be calculated in §2.4[ 

In the following we shall omit the index tot for semplicity, while we will 
always explicitely write down the index when considering a specific species (for 
instance ge stands for electron energy DOFs). 

Time - temperature relationship 

From the definition if = - we get: 

t = -—. (2.45) 

JO Ha 

By using the equations fl2.12p . neglecting the curvature term in the radiation 
dominated era and recalling equation (12.401) we get: 



H ^ J^^g^^T' ^ IMg'/'^, (2.46) 



where we used G ~ Mp^. From (12.461) we can derive the relationship between 
the time t and the background (photon) temperature T, using the entropy 
conservation and assuming both g and q approximately constant. This way 
we find 

-1/2 Mp _ f T 



t ^ 0.301 g~'^' ^ ^ t^ [j^] (sec) (2.47) 
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2.3.2 Thermal equilibrium criterion 

The correct way to evolve particle distributions is to integrate the Boltzmann 
equation; nevetheless a rough criterion to evaluate whether a particle species is 
in thermal equilibrium or not can be found considering that a species which is 
not in equilibrium at time t will never thermalize if every particle has in average 
less than one interaction from t to cxd. 

Let us consider two interacting particles species, one of which will be treated 
as a target. The interaction rate is defined as: 

r = n{va) (2.48) 

where n is the number density of terget particles, v is the relative velocity and a 
is the interaction cross section. If the particles are in thermal equilibrium with 
the photon bath, we can assume that F is some power of the photon temperature 
T; therefore, using equation fl2.47p we have: 

r oc T'^ oc t-"/2 (2.49) 

From equation (12.461) we can see that the Hubble constant H oc t^^. Therefore 
the number of interactions per particle for time going from t to oo is: 



T{t') dt' oc T{t'T dt' 



oc 



dt' oc ^ _L (2.50) 

H 



Therefore the commonly used rule: 

r(T) > H(T) =^ implies thermal equilibrium (2-51) 
from which we can define the decoupling temperature Tjj as: 

"Mr' <-^) 

Note that T{T) < H{T) does not imply departure from thermal equilibrium be- 
cause a non-interacting species once in equilibrium will ever keep an equilibrium 
distribution. 



2.4 Neutrino decoupling and e annihilation 

In the early universe neutrinos are kept in equilibrium via the reactions z/z/ 
e~^e~, PC vc etc. Let us assume that all the involved particles are ultra- 
relativistic. The cross sections - see §B.1I - are of order: 

a ~ Gls ~ GlT^ (2.53) 
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where Gp is the Fermi constant Gp — 1.1664-10 ^GeV ^ and we used s oc as 
in §A.5[ Since n c:^ and c = 1 the interaction rate defined in (I2.48P becomes: 

r ~ GlT^ (2.54) 

Using (12.461) and (I2.54p in the definition of the decouphng temperature given in 
(I2.52P we can see that the neutrinos decouple at the temperature Tjj given by: 



F-^D 



H{Tn) Tl/M, 



^ IMeV (2.55) 



Therefore at T ^ IMeV neutrinos are in thermal equilibrium with the plasma 
and their temperature is T^, = T, while after decoupling their temperature scales 
as a~^. 

Shortly after the u decoupling, the annihilate because T becomes smaller 
than 2me, which is the threshold for the reaction 7 ^ e~^e~ . Thus electrons and 
positrons transfer their entropy to photons, which become hotter than neutrinos. 

A more quantitative under standment of these temperature relationships will 
be necessary in chapters [5] and [6] and can be found below. 

Neutrino temperature after decoupling 

After the z/s decoupling, the i/s' and the remenant ultrarelativistic particles' 
entropies S are separately conserved: 

Siy(T) = Sy{T) ■ = const = ci 

Sje±{T) = s^e±{T) ■ = const = C2 (2.56) 

The ratio of these entropies is also a constant, which we will call C3; in the ratio 
the acceleration factor simplyfies, giving 

'"^'^"^ - const = cs (2.57) 



S^e±{T) 

When the annihilation takes place, photons and neutrinos are ultrarelativis- 
tic; therefore their entopy densities can be worked out from equations (12.381) 
and ([232]): 

27r2 7 n 7 . 
27r2 4 

s,{T) = —2T' = -asT' (2.58) 
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where T is the photon temperature and ob is the radiation constant: 

as ■■ 



— (}i = c=l units) 
15 



(SI units) 



(2.59) 



15(|(c)3 

Electrons and positrons are non relativistic but in equihbrium with the photons; 
their entropy density can be therefore worked out using equations (12.341) . (12.351) 
and fl236|) : 



Se±(T) 



n^T \ 3 



m 



2^3 



dE + 

•e exp T +1 Jnie 



oo J{E^ - m?)E^ 



exp T +i 



dE (2.60) 



where we assumed /ie± = 0. Equation (12.571) finally becomes, using the entropy 
densities in f l238|) and ^M)- 



2T^ + Se±{T) 



15 



ca 3 



7 ^3 1 27r2 
2^^ '^"^[l^ 
where the function / is 

and the new variables x and y are defined as 



1 3 



3x2 _^ 4^2 



me 
T 



= ^ 



(2.61) 



(2.62) 



(2.63) 



Using the same formalism it is straightforward to write the electron-positron 
energy density as: 



2 fo^ ^{E^-m^)E^ ,30 f^^/^+¥}y^ 
Pe± = — I — — : — dE = aBT — ^== dy (2.64) 



2 / -E-M 
71 Jrue exp T +1 



vr^Jo expV^ 



+1 



and therefore the total energy density as a function of x and T: 



p(T, x) 



'3,2 _^ ^2)^2 



(2.65) 



exp V ' w +1 

The first term is the photon's energy, the second one is the neutrinos' and the 
third is the electron-positron's. The neutrino term is weighted by the number 
of degrees of freedom ^ and the temperature factor. 



Equation (I2.6ip totally determines the neutrino temperature as a func- 
tion of the photon temperature T once the initial conditions are specified by 
mean of the costant C3. 

It is worth to stress that the status of the annihilation process is totally con- 
tained in the / function as the ratio x = Its extreme cases are: 
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• X <^ 1 =^ T ^ me, that is ultra-relativistic electron. In this hmit the 
reactions that convert in photons and viceversa are in equihbrium and 
the annihilation process is not begun yet. In this limit / becomes: 

+ — ^ '-dy^- (2.66) 

• X ^ 1 =^ T <^ mg. The leads the integral to giving 

[/(x » 1)]=^ ~ 1 (2.67) 

This limit describes the situation we have after the end of the annihilation 
process. 

Working out the constant c^: the instantaneous decoupling approxi- 
mation 

Let us assume that the neutrino decoupling is an instantaneous process taking 
place when photons have the temperature T = Tn^, before the decoupling, 
neutrinos and photons have the same temperature and therefore 

_ lqe{T) + q^ _ lqei{TD,) + 2 

^3 - - — ,3 — y^-^^) 



^qy\T ) 
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The value of the constant C3 is affected by the difference of T^j, and the 
annihilation temperature Tann- If these temperatures are close indeed and the 

annihilation begins before the v decoupling is completed, the transfer a 
fraction of their entropy to neutrinos, raising their temperature. 

The asymptotic values of C3 are: 

1-4 + 2 22 

ca = ^ 7 ^ = — = 1.048 Tdu > Tann 
g - 6 21 

C3 = = |- = 0.38 Tdu < Tann (2.69) 

g - 6 21 

Tipical values for T^y are 1 4- 3 MeV\ numerical calculations show that in this 
range C3 = H within 2% of error. Therefore in chapters O and [H] we will just use 

C3 = ^ 
^3 21- 

Using this value for C3 we also get the standard asymptotic ratio of the 
neutrino and the photon temperature: 

T, = ' T (2.70) 
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which is valid when T <^Tc 



ann- 
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2.5 Primordial nucleosynthesis (BBN) 

The Big Bang Nucleosynthesis (BBN) is the earhest test of the cosmological 
and of the particle interaction models. Energy considerations suggest that light 
nuclei could be formed when the temperature of the universe is in the range 
1 — SOMeV; nevertheless, we will demonstrate below that, because of the high 
entropy of the universe, the actual synthesis takes place at a much lower tem- 
perature, T/v ~ O.lMeV. 

Let us start introducing the mass fraction of a certain nuclear species, which 
is defined as 

Xa^^ , E^^ = l (2-71) 

where A is the atomic number, stands for nucleons and ua is the number 
density- see §2.31 and §A.3[ Xa can be recasted as 

2 A~l 

X^ = g^Ah''^Ct'n'^(—) ' v^~'X^Xt'exp^ (2.72) 

where we neglected the difference in mass of proton (p) and neutron (n). In 
equation fl2.72p Ba is the binding energy of the nuclear species A, defined as: 

Ba = Zrrip + [A — Z)mn — rriA (2.73) 

This quantity has values which vary from 2.22MeV (for ^H) to 92.2Mel^ (for 
^^C), corresponding to a binding energy per nucleon of order of 1 to 8Me\^; 
finally rj is the barion to photon ratio in the universe and is proportional to s~^: 

riB vr"^ nBq{T) ^nBq{T) 

r] = — = ~ 1.8 ^ — - 

rij 45^3 s s 

riB = rib-n^-riN,^,^^ (2.74) 

The main elements produced during BBN are D, ^He, '^He, Li with a predom- 
inance of ^He, having a mass fraction of about 24% at the end of BBN. Among 
the light elements, D and ^He play a crucial role because there are apparently 
no astrophysical processes that can account for their observed abundances; it is 
therefore nowadays possible to measure to high precision their primordial abun- 
dances and, by comparing the measured and the theoretical values, it is possible 
to test the cosmological model and also to work out bounds on MCP - see §4.51 
and chapter O 

We can infer from equation (I2.72p that, although the binding energies per 
nucleon are 1 -h 8MeV^, the equilibrium abundance of nuclear species do not 
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become of order unity until a temperature of order O.SMeV because the high 
entropy of the universe leads to low values of r] - see equation (12.741) . A rough 
estimate of when a species becomes thermodinamically favored can be worked 
out assuming X„ ~ Xp ~ 1 and solving for ~ 1, from which follows 

- (2,75) 



— \n{rj) + 1.51n(mAr/T) 

This temperature is much lower than the binding energy per nucleon Ba/A: for 
instance 

= = 0.07MeV — = l.llMeV 

2 

T3He = O.llMeV = 2.57MeV 

T4He = 0.28MeV = 7.07MeV 

4 

Ti2c = 0.25MeV = 7.68MeV (2.76) 

Therefore the high universe entropy favour free nucleons for T < T^; since 
these temperatures are much lower than the nuclear statistic equilibrium ones, 
BBN can not begin when the species go out of equilibrium, but has to wait till 
T ~ T4. Moreover, the reactions producing helium and heavier elements are all 
based on D: 

D + D ^He + n 
^He + D ^He + p (2.77) 



D + D ^H + p 
^H + D ^He + n (2.7J 



D + D *He + j (2.79) 

and thus the reactions can not be fast enough to produce an equilibrium abun- 
dance of ^He for T < O.lMeV. 

Protons and neutrons are kept in chemical equilibrium by mean of the weak 
interactions: 

n ^ p + e~ + 
V + n ^ p + e~ 
n + e+ ^ p + i^e (2.80) 



§2.51 Primordial nucleosynthesis (BBN) 



33 



These processes have approximately the rate Tw — GpT^ and go out of equi- 
hbrium at the "freeze-out" temperature of weak interactions, when T(Tw) — 
H{Tw), that is: 

Tw ^ (0.7 - 0.8)MeV {tw ~ Is) (2.81) 

The neutron to proton density ratio, that is ^ = is given, when chemical 
equilibrium holds, by 



n Tin 

— = — ^ exp 

p n. 



-Q 



T 



op 

where Q is the neutron-proton mass difference: 



(2.82) 



Q = m„-mp~ 1.293 Mel/ (2.83) 

For T < Tw the weak reaction rate Tw drops below the Hubble expansion rate 
H{T) ^ b.bT"^ /Mpi, the neutron abundance freezes out at the equilibrium value 
Xn{Tw) and it then evolves following the neutron decay exponential law: 

X„(t) = X^{Tw) exp(-t/r) (2.84) 

where r = (886.7 ± 1.9)s is the neutron lifetime. At this time the neutron to 
proton density ratio is 

-{Tw) ^ I (2.85) 
p o 

This quantity further decreases to 

-(Tn) ^ ^ (2.86) 
p 7 

due to the occasional weak interactions, mainly neutron decays. At tempera- 
tures T < Tjv ~ O.lMeV, the process p + ri'^D + 'j is faster than the Universe 
expansion, and free nucleons and deuterium nuclei are in chemical equilibrium. 
Also, reactions are fast enough to produce an equilibrium abundance of '^He, 
which can be worked out assuming that all neutrons end up in '^He. In conclu- 
sion, the primordial ^He mass fraction is 

_ 4{nj2) _ 2H(T^) _ 2exp(-Wr) , 
' " n„ + - 1 + n(T;v) " 1 + exp{Q/Tw) " ' ' ^ 

The ^He production is influenced by the number of degrees of freedom g, which 
enters H, by the neutron half life r and by t]. Comparing the theoretical calcula- 
tion of the primordial '^He abundance, today accurate to within ±0.4% [78] with 
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the experimental observations - concordant within ±4% [80] - it is possible to 
bound these three parameters. In particular, the number of degrees of freedom 
can not be higher than a certain value, otherwise too much ^He is produced. 
It is common to re-parametrize the DOFs number in terms of extra-neutrinos, 
defined as 




(2.88) 



Bounds on AA^^, have been calculated by several authors - for instance Lisi et 
al. [77], who obtained 



AN^ = ± 1 



(95% C.L.) 



(2.89) 



Chapter 3 



Electric charge quantization in 
theoretical physics 



The puzzle concerning the electric charge quantization has been studied since 
the beginning of the 20^^ century and is nowadays still challenging. 

There are several experimental and theoretical hints suggesting that the 
electric charges of the particles entering the SM are quantized. For instance 
measures of the neutrality of matter from binary pulsars [23] and magnetic 
effects [24] set respectively 



where and Qe are the proton and electron electric charges, while the neutron 
one Qn is [25] 



Other interesting measures concern the muon-electron charge ratio anomaly 
^ + 1 = (1.1 ± 2.1) X 10-9 [26] and the neutrino charge < 10"^^ [27]. 
We defined in §1.51 the SM electric charge in the SM scenario as 



This quantity is not quantized in the classical field theory since the hypercharge 
is not. Nevertheless, we will show in §3.1l that imposing the SM renormalizability 
at quantum level we will get for each generation four relations on the five present 
hypercharges. This topic can be found in many books about quantum field 
theory, such as [28-30]. 

Also several models beyond the standard one provide mechanisms to quan- 
tize the electric charge. Some of them achieve quantization adding new particles 




(3.1) 



= (-0.4 ± 1.1) X 10~2ig 



(3.2) 



Qem ^3 ~l~ 



(3.3) 
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or imposing features to the existing ones (such as the neutrino being a Majo- 
rana particle in [35,36]) or using gauge groups different from the SM one (for 
instance [5W(3) x SU{2) x U{1)]l x [SU{3) x 5W(2) x U{1)]r in [34]). 

Electric charge quantization also arises as a natural consequence in Grand 
Unification Theories and in Dirac quantistic theory of magnetic monopoles, 
analyzed in §3.21 and 13.31 These topics can be also found in textbooks, see for 
instance [29,30]. 

We conclude this brief review stressing that, in spite of these constraints, 
there are mechanisms, like the photon kinetic mixing introduced in §4.H which 
are able to add particles with unquantized charge to SM extensions gauged by 
two and to some gauge theories without spoil their renormalizability or its 
other features. 

3.1 Anomalies and their cancellation 

An anomaly is the failure of a classical symmetry of the lagrangian C to survive 
the process of quantization and regularization. Indeed, if we have a classic 
symmetry, the transformation (p ^ (p + 6(j) will leave the action S{(j)) invariant, 
while, if we have a quantum symmetry, the same transformation will leave the 
path integral / Dcpe^^^'^^ invariant, where Dcp is the measure. Therefore it looks 
reasonable that some classic symmetries may be not valid in quantum theories. 

An important example of anomaly concerns the chiral symmetry of theories 
with massless fermions. Before entering this topic, let us define the vector, 
vector axial and pseudoscalar currents, which are respectively: 

V^'{x) = V^(x)7^V^(x) (3.4) 
A^(a;) = V^(x)7'^7^^(x) (3.5) 

P{x) =^{x)-f^ip{x) (3.6) 

We will see that the currents (13. 4p and (13.51) are conserved in the classic theory, 
but it is impossible to preserve both these conservations in the quantum one. 

3.1.1 Chiral symmetry for free Dirac spinors 

Dirac spinors ip and their adjoints ip = ?/^^7° obey the Dirac equations: 



— m)ilj = + m) = 



(3.7) 



§3. Jl Anomalies and their cancellation 



37 



which can be derived from the Dirac lagrangian 

jO-Dirac = - m)ilj (3.8) 

The Dirac lagrangian fl3.8p is manifestly invariant under the phase transforma- 
tion: 

^(x) ^ e'^'tpix) (3.9) 

Since this continuous symmetry of £, there must exist for the Noether theorem 
- see §1.11 - the conserved current defined in equation (11. 6p . For the symmetry 
we are analyzing the Noether current is the same than in equation (13. 4p and 
thus 

d^Vix) = (3.10) 

The conservation of this current can be also derived using directly the Dirac 
equations (13. 7p . Let us consider now the axial vector current (13. 5p : using the 
Dirac equations it is straightforward that: 

d.Af'ix) = 2imP(x) lim d.j^^x) = (3.11) 

m — >0 

Thus, when fermions are massless, the theory has a second Noether current 
corresponding to the chiral transformation: 

V'(x) ^ e^^^'V'la;) (3.12) 

The chiral symmetry pertains to the only derivative term in the lagrangian (13. 8p : 
this is why the axial vector current is not conserved for massive fermions. 

When m = the currents in (13.41) and (13. 5p can be used to define the electric 
current densities of left- and right-handed particles, which are both conserved: 

Mx) = ^{x)Y i^^r^] i^ix) = I [V'ix) - A^{x)] ^ ^^r^L (3.13) 



In conclusion, the massless Dirac lagrangian has a symmetry associated with 
the separate number conservation of left- and right-handed fermions, generated 
by the axial vector current j^^. 

The conservation laws presented above are also valid in the quantum theory 
of free Dirac spinors. 
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3.1.2 Chiral anomaly in abelian gauge theories 

In gauge theories interactions are present - see §1.21 When quantizing the theory, 
the functions we have in the classic model are replaced by quantum operators; 
if we assume these operators as a simple generalization of the classic functions, 
they may be not well defined and thus behave in a different way with respect 
to the classic ones. In this section we will see that one of the arising effects is 
the spoiling of the axial vector current conservation law. 

Let us consider now QED, which is an abelian gauge theory - see §1.2.11 -: 
its lagrangian is0: 

Cqed = i^{ip-m)ij--F^,F^''' (3.15) 



D>' = d''- ieA'' F"'' = d^A" - d^A^" (3.16) 

The lagrangian in ( 13.15P does not change under the local gauge transformations: 

^ ^ exp^^°(^) A^^A^- ied^a{x) (3.17) 

As it happens for free spinors, if m = the lagrangian is also invariant under 
the local axial symmetry: 

7/- ^ exp^"°(^)^' ^/^ A^^ A^-ied^a{x)-i^ (3.18) 

which is broken by the mass term. Therefore we would naively expect d^j^^ = 
in massless gauge theories. 

Nevetheless, performing the calculations more carefully, we can see that the 
actual picture is more complicated in the quantum scenario. To get an intuitive 
idea, recall the equal-time anticommutation relationship for Dirac particles in 
second quantization: 

{U^),4{y)} = 6,,6^'\^-y) (3.19) 

When X = y this anticommutator diverges; thus we may wonder if the quantity 
ijj{x.)ilj{x.) is physically meaningless or not. Because this divergence, also the 
classic currents defined in (13. 4p . (13. 5p and (13.61) . which contain the product 
'0(x)?/'(x), may get problems when generalized to quantum operators. 

What comes out indeed is the impossibility to preserve both vector and 
axial vector current conservation. Both these currents finally depend on an 
arbitrary parameter, which is chosen such that the gauge vector invariance is 

-•^In this section we use the same conventions of [28]; several books, such as [31] define F'^'^ 
and other quantities with opposite sign. 
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preserved and the theory of electromagnetic interactions with massless photon 
is safe. But this also implies that the derivative of the axial vector current gets 
non-zero value in the massless limit§|: 

limd^A^ = --^e-'^M-F^^F^, (3.20) 

Using the current densities of left- and right-handed particles introduced in 
eq. (13.13p and (13.14p we can rewrite the anomaly for massless particles as: 

1 p2 

In conclusion, the massless lagrangian C is invariant under the chiral trans- 
formation in eqs. fl3.18l) . but the associated current is non conserved. This 
property of the quantum gauge theories is called Adler-Bell-Jackiw (or ABJ) 
anomaly [40,41]. Its validity to all orders in QED perturbation theory has been 
proved by Adler and Bardeen [42]. 



3.1.3 Chiral anomaly in non-abelian gauge theories 

Non-abelian gauge theories are also anomalous as we can easily demonstrate 
once that equation fl3.20l) is assumed. Let us consider a gauge theory with a 
group having a set of generators and the commutation relation between them 
in the form: 

[t'',t^]=ir^H'' (3.22) 

where f"''"^ is an antisymmetric set of numbers called structure constants. In 
this picture a theory containing massless Dirac fermions has the lagrangian: 

The same mathematical procedure leading to equation (13.201) can be applied to 
the gauge symmetry current: 

/'^=?7^(^^)t> (3.24) 
2^ap(ii> jg j^Yie totally antisymmetric tensor defined such that e"^^'^ = — eoi23 = +1- 



d. - igAT 



r 



(3.23) 
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Assuming massless particles we finally get [43]: 

{p,u,b-k,X,c\d,jnO) = ^^e^'^'p.kf, (itr[tnt^t^}]) (3.25) 

Of course, if right-handed fermions enter the triangle loop there will be analo- 
gous terms but with opposite sign. It can be shown that the total anomaly is 




Figure 3.1: Triangle diagrams contributing to the axial vector 
anomaly. 



proportional to the one arising at triangle level, where the anticommutator in 
the trace is easily explained if we look at the diagrams in Figure 13.11 the total 
triangular anomaly is worked out summing the two diagrams having {p, u) and 
{k, A) interchanged and this leads to the anticommutator in the trace. 

The factor | in front of the trace is a kind of weight which can be understood 
requiring that if we take = t^ = = 1, we should get QED and therefore the 
term between round brackets should be 1. Indeed: 

^tr[l{l,l}] = l (3.26) 

Physically consistent gauge theories must be anomaly free. From equation (13.250 
we can see that this condition is achieved if the generators satisfy 

tr[e{t^t=}] = (3.27) 

But there are also anomaly-free theories containing anomalous diagrams which 
cancel each other. In §3.1.5l we will see that one of these theories is the standard 
model of elementary particles. 

3.1.4 Chiral symmetries in QCD and the tt^ decay 

Let us consider QCD with the only u and d quarks; its lagrangian is: 



C = uipu + dipd — rriuUU — madd 



(3.28) 
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In the following we will neglet the quark masses. Besides the U{l)v x W(1)a sym- 
metry, where V and A stands respectively for vector and axial, this lagrangian 
is symmetric under the unitary transformations: 

because there are no couplings between left- and right-handed quarks. Let 
Q = Ql + Qb. denote the quark doublet, with chiral components: 

The theory is symmetric under the group SU{2)l x SU{2)r x U{1)l x U{1)r, 
which corresponds to the conserved currents: 

fL = QlI't'^Ql fR=QRYr''QR (3.31) 

with t"" = ^ generators of SU{2). Equivalently we can say that the theory 
is symmetric under the group SU{2)y x SU{2)a x U{1)v x U{1)a] sums and 
differences of the currents in (13.311) indeed lead to: 

Af" = Q-i^-i^Q A^" = Q7Vr"Q (3.32) 

Classically the four currents in fl3.32p are conserved; nevertheless, when passing 
to the quantum theory, we expect anomalies to arise for axial currents - see 

An important application of the ABJ anomaly in QCD is in the derivation of 
the theorem for the vr'' 27 decay in two photons - see figure 13.21 Applying the 
formalism we developed in §3.1.3[ it is straightforward to see that this process is 
described by a triangle diagram having the axial isospin current anomaly given 
by dS^OD, multiplied by 

^tr[r^{Q,Q}]=tr[r^Q2] (3.33) 
where Q is the matrix of quark electric charges: 



(3.34) 
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Figure 3.2: Triangle diagrams contributing to the 7r° decay in two 
photons. 



The trace runs over flavours and colors; the involved matrices do not depend 
on color, therefore we just get from it a factor 3. Since tr[r^Q^] = tr[r^Q^] = 0, 
the axial isospin current reads: 



If there had not been anomalies, the amplitude for the vr" decay would have 
been 0, that is, not compatible with experiments. Moreover, the decay rate F 
is related to the coefficient of the ABJ anomaly, which can be this way experi- 
mentally tested. The predicted rate 



where /tt = 92.4 MeV is in very good agreement with the measured value 
F = 7.7 ±0.6 eV [86]. 

3.1.5 Anomaly cancellation in the standard model 

We have demonstrated above that the presence of anomalies is not forbidden, 
on the contrary these diagrams can sometimes be useful. But they also spoil the 
gauge invariance of the theory; therefore the anomalous terms must cancel each 
other in consistent gauge theories. Some models anomaly free can be found 
in [44,45], while the general features of a gauge theory free of anomalies are 
analyzed in [46]. 

We will show that the electro-weak model first proposed by Glashow-Weinberg- 
Salam [1-3] is not anomaly free, while the whole SM SU{3) x SU{2) x U{1), 
in which quarks are weighted by the color factor 3, does. This is a remarkable 
proof of internal consistency of the SM. 



167r2 



(3.35) 




(3.36) 



§3. Jl Anomalies and their cancellation 



43 



Electroweak anomalies 

Let us consider now the only electroweak sector, that is the gauge group SU (2) x 
To evaluate its anomalies it is easiest to work in the basis of the gauge 
bosons before the mixing leading to the photon and definition - see equa- 
tion fll.SSp . We therefore have four generators, three from SU{2), which are 
proportional to the Pauli matrices at, and one from the weak hypercharge Y: 




= 1-3 



Y = ^ (3.37) 

The Pauli matrices obey to 

{cTj, aj} = 2Sij tr[ai] = (TiCTj = ie^kcrk + 5ijl2x2 (3.38) 

where eijk is the antisymmetric Levi-Civita symbol. Let us analyze now all the 
possible combinations. Equations (13.381) imply that 

tr[ri{rj,rk}] = ^5jk tr[ai] = (3.39) 

Considering then that every member of a SU{2) multiplet has the same hyper- 
charge, we get: 

tr[ri{Y,Y}] octr[ri] = (3.40) 

In general, any anomalies containing one only SU(2) boson are proportional to 
tr(ri) ans therefore vanish. There are therefore only three terms which do not 
automatically vanish. The first concerns the left-handed particle hypercharges: 

tr[Y{rj,rk}] = ^5jk tr[Y]L = ^^jk ( E Y + E Y ) (3.41) 

yieptons quarks J ^ 

From which follows the condition to have an anomaly-free model 

3y,, + F/, = (3.42) 

where the factor 3 arises summing on the three quark colors. The second lead 
to the same condition by using the third equation in fl3.38p in the trace 

tr[ri{Y, Tj}] = 2tr[rirjY] ex tr[(ieijkak + 5ijl2x2)Y] cx (5jktr[Y] 

=^ tr[ri{Y,rj}] = ^5jktr[Y]L (3.43) 
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The third term is tr[YYY] and lead to a different condition: 

tr[Y3]=0 =^ ^y3 = (3.44) 

which becomes, if written exphcitely for the first particle generation: 

6Yl + 2Y^l - (3F4 + 3Yl + vQ = (3.45) 

In conclusion, the two equations (13.421) and (13.451) are enough to ensure the 
anomaly cancellation in the electro-weak sector of the SM. 

Color anomalies 

Let us consider now also the color interactions; the SU{3) generators are the 
eight Gell-Mann matrices A**, having the properties: 

tr(Aa) = 



tr(AaAb) = 2S, 



ab 



4 

{Aa, Afe} = - 5a{) l3x3 + 2rf„fec Ac (3.46) 



where dahc is the totally symmetric symbol defined as: 

1 



<^118 — ^^228 — "^338 



(i448 — (isss — (i668 — '^778 — 7:^^888 



2 2^3 

C?i46 = (ii57 = — (i247 = C^256 = = (^355 = ~C^366 = ~C^377 = ^ (3-47) 

QCD is a left-right symmetric theory, such as QED, and therefore its anomalies 
cancel each other. Also, traces containing only one SU{2) or SU{3) bosons are 
proportional to tr{Ti) and tr(ta) and therefore vanish. The only terms which 
survives is 

tr[Y{Aj, Ak}] = tr[Y] = J^h^Y ^'^jkY^ (3.48) 

where the sum runs over left- and right-handed quarks and / is 1 for left-handed 
and for right-handed quarks. To get anomaly cancellation therefore the model 
must satisfy 

2Y^^ - r., - Irf, = (3.49) 
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Gravitational anomaly 

There is also a gravitational anomaly with one U{1) which is proportional to 
tr[Y] and leads therefore to the condition 

tr[Y] = ^y = 



=^ 2Yi^ - Yi^ + QY,^ - - 3F,^ = (3.50) 

summed over left- and right-handed leptons and quarks. 

In conclusion, the five hypercharges we have for any particle generations 
must satisfy four conditions to lead to the anomaly cancellation in the SM: 

3y,, + i^, = 

6F/^ + 2y,J - (3F4 + 3it + = 



2Y —Y —Yri 



R 



2y,, - r,, + er,, - 3K, - 3r,, = o (3.51) 

These equations can be recasted in terms of the electric charge using the electric 
charge definition (13.31) . from which follows 

Y = 2{Q- rs) (3.52) 

and the is tracelessness, from which tr[Y] = 2tr[Q]. 



3.2 Grand unification and SU{5) 

Grand Unification (or Unified) Theories (GUTs) are based on the intriguing idea 
that at extremely high energies, of order 10^^ GeV, the three forces involved in 
the SM are unified in a single group with one coupling constant. At low energies 
this group is broken in the familiar SU{3) x SU{2) x U{1). 

Leptons and quarks are tipically put together in GUT multiplets; as a conse- 
quence, quantization rules for the electric charge naturally arise in these theories. 
Some models trying to evade this feature can be found in literature - see for 
example [33] - but they need the addition of several new particles. 

It has been shown by Georgi and Glashow in 1974 that the simplest grand 
unified model including all the SM features is SU{5) [39]. Indeed the unified 
group should: 
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• Be at least of rank 4 because the SM has four mutually commuting gen- 
erators, one from W(l), one from SU{2) and two from SU{3). 

• Have complex representations, because in the SM fermions are not equiv- 
alent to their complex conjugates, that is they tranform differently. 

• Take in account of the existence of both integer and fractional charges. 

Besides having all the features written above, SU{5) has the anomaly-free repre- 
sentation 5-1-10 , in which the observed quarks and leptons fit neatly getting the 
correct quantum numbers. In the following this representation and its particle 
content are analyzed and the electric charge quantization rule is shown. 

3.2.1 Representations and particle content - one family 
case 

Every lepton family entering the SM contains 15 chiral modes; let us consider 
for semplicity one lepton family e and its chiral modes, which are listed below 
with their transformation properties under SU{3) x SU{2) x U{1): 

{ue,e-)L : (1,2,-1) 

{Ua , dc,)L ■ ' ^ ' ^ 

(Ol : (3,1,-^ 
{d^L : (3,1,^) (3.53) 

The W(l) transformation properties will be omitted in the following. The su- 
perscript c indicates the charge conjugate field: 

= = C*^^ (3.54) 

SU{5) does not have a 15- dimensional representation, thus the fermion content 
is split into the sum of a 5-dimensional and a 10-dimensional representation. 
The 5 is assumed to be right handed and the 10 to be left handed, so that 
the anomaly-free combination 5 -|- 10 is left handed and contains all the modes 
hsted in ( l333|l . 
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The fermion content of this representation is therefore: 
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(3.55) 



(3,1)©(1,2) iV^g: 
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(3.57) 



The gauge mesons transform according to the adjoint representation of SL({5), 
which has dimension 5^ — 1 = 24 and can be decomposed according to the 
SU{3) X SU{2) quantum numbers as: 

24 = (8, 1) © (1, 3) © (1, 1) © (3, 2) © (3, 2) (3.58) 

The gauge bosons are identified as follows: 

• (8, 1) are the eight colored bosons of SlA{?>)c] 

• (1, 3) © (1, 1) are the four SU{2) x U{1) bosons, that are 7, and Z°; 

• (3, 2) © (3, 2) are new superheavy gauge bosons which couple the quarks 
to the leptons and mediate the proton decay. They are usually denoted 
as the X, Y bosons. 



3.2.2 Three families generalization 

In the one family approximation gauge and mass eigenstates coincide. To gen- 
eralize to three families we have to replace the fermion fields with gauge eigen- 
states, which are vectors in the three-dimensional space of the family index 
A = e , 11 , t: 

e ^ = ^AB^B eB = {e,fi,T) 
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u-^Pa = U\b Pb Pb = (u, c, t) 
n'j^ = v\b ub ub = {d, s, b) 



- P'l 



U%Pb PB = {u'',c\t 



C „C 4.C\ 



d'-^n'X = YXbTi^ nB = id',s',b'') (3.59) 

The basis has been chosed such that lepton-gauge and mass eigenstate ca co- 
incide. For neutrinos we should generally have instead of 6ab, but if we 
assume that neutrinos are massless, any linear combination of the degenerate 
fields can be taken to be their mass eigenstates. 

If we wish to work only with left handed fermions, the transformations for 
left- and right- handed particles should be different. This is why we wrote down 
exphcitly the last two equations in fl3.59p where 

Uab + U%B Vab + VXb (3.60) 



3.2.3 Charge quantization 

The charge quantization arises in the SU{5) scenario because this group is 
simple and non abelian. Q is an additive quantum number and has therefore to 
be some linear combination of the four diagonal SU{5) generators. 

Among the 24 SU{5) generators, eleven are just the generalization in five 
dimension of the SU{3) and SU{2) generators, namely the Gell-Mann and Pauli 
matrices: 

L'^ = I ° ^; a=l,2,...,8 

= 3 ); 6 = 9,10,11 2 = 1,2,3 (3.61) 

from which we get three diagonal generators: 



/ 1 











\ 





-1 









































lo 











0^ 



§3.21 Grand uniGcation and SU (5) 



49 



2-^3 
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(3.62) 

1 

V -1 y 

The fourth diagonal SU{5) generator is a traceless matrix which commutes with 
the others; since the identity matrix commutes both with SU{3) and SU{2) 
generators, this fourth diagonal generator is split in two diagonal blocks, both 
proportional to the identity in their dimension, that is 



12 



a 



Isxs 
/3 





l2x2 



(3.63) 



From the traceless and the normalization condition TrL'^V' 

T 



2(5"^ we get: 



3a + 2/5 = 
3^2 + 2/52 = 2 



and therefore 
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(3.64) 



(3.65) 



As stated before, Q must be a linear combination of the diagonal generators in 
f l3:62|) and fl335D . that is 



Q = aL^ + [3L^ + 7L^^ + 5L^'^ 
Let us now apply Q to the fermion representations in ( 13. 55^ : 
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Therefore we can deduce that the electric charge is given by: 
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From the traceless of Q we get, for ric = 3 colors 











(3.67) 



(3.68) 



ric- Qd + Qe + Qu = ^ Qd = ~^Qe 



(3.69) 



Thus the electric charge quantization arises in SU{'b) as a consequence of having 
three quark colors. 



3.3 Magnetic monopoles 

The theory of magnetic monopoles was first formulated by Dirac in 1931 [47]. 
As it is shown in the following, it is consistent with quantum theory only if the 
electric charge is quantized. 

3.3.1 Dirac classic theory of magnetic monopoles 

The original Dirac's theory is a straighforward generalization of the classical 
electromagnetism, which is modified to include magnetic sources. Classical 
electromagnetism is based on the Maxwell's equations: 

V-E = p V-B = 

rfVK riV 

VxE = VxB=j + - (3.70) 
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the solutions of which, E and B, can be parametrized in terms of the vector 
and scalar potential A and as: 

OA 

E = -V0-— B = VxA (3.71) 

It is worth to stress that we can write the second equation in (13.711) only in 
consequence of the Maxwell's equation stating V • B = 0, because for any- 
arbitrary vector V we have V ■ (V x V = 0). 
Maxwell's equations can be written in the covariant form: 

duF^'' = (3.72) 
in terms of the electromagnetic field tensor F^^: 

F^, = d^A,-d,A^ A^ = {<P,A) (3.73) 

and its dual F^'^: 



2 



'Fp^ (3.74) 



The only source in f l3.72p is the electric current j'^ = (p, j); in its absence, that 
is in vacuum. Maxwell's equations are symmetric under the duality transforma- 
tion: 

This symmetry is broken by the current j^, but it can be restored introducing 
the magnetic current k^^ = (ex, k) in eqs. (13.721) . so that the covariant Maxwell's 
equations become: 



d^F^"" = -J 



d^F^"" = (3.76) 

and the duality symmetry holds if the substitution in (I3.75P are made together 
with: 

f ^k" F ^ -j" (3.77) 

In classic electromagnetism the current is produced by electrically charged 
particles; similarly the magnetic current k^^ can exist only if there are magneti- 
cally charged particles called magnetic monopoles. 
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3.3.2 Quantum theory of magnetic monopoles and elec- 
tric charge quantization 

To work out the quantum behaviour of Dirac's magnetic monopoles, we recall 
equations fl3.76p and their solutions for point electric charge/magnetic monopole 
fixed at the origin: 

V • E = 47rg S^{r) ^ = 

V • B = 47r^ 6^{r) ^ = 9^ (3-78) 

Let us now consider the surface integral for B over any closed surface S con- 
taining the monopole: 

jB-dS = ATi-g (3.79) 

This integral is not zero because V • B ^ 0; as a consequence, we have for 
arbitrary x, B(x) 7^ V x A(x) as we stressed below eqs. f l3.7ip . 

A trick to evade this condition can be found considering the field generated 
by a long and thin solenoid placed along the 2;-axis with its positive pole placed 
at the origin. The solenoid's field is: 

"Qsoi = + ge{-z)5{x)5{y)i (3.80) 

The solenoid's field has no magnetic sources, therefore we can write: 

= V X A (3.81) 

Comparing fl3.78p and fl3.80p we see that 

^monopole = B = B^^, - gO{- z)5{x)5{y)z 

= V X A - ge{-z)5{x)S{y)z (3.82) 

that is, it is possible to define A such that B = V x A everywhere except a line 
going from the origin to infinity, which is called Dirac string. 

The quantization of electric charge follows requiring that the solution of the 
system with the monopole is mono-valued when integrating around a loop. 

The solution of fl3.78p is a plane wave: 

= C-exp^(P-'-^*) (3.83) 

with the substitution: 



p^p-qA ^ V ^ "0 ■ exp *'?(^-'') 



(3.84) 
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which leads to an additional phase factor. 

The wave function is mono-valued when performing an integration around 
a loop only if its phase factor has the value exp*^™ = 1, which is achieved if: 

27m = q(j^A-d\ (3.85) 
Using flXT^ and flXS^ we finally have: 

q j A- d\ = q j B ■ dS 2nn = q4:ng (3.86) 

The quantization condition induced by the presence of the Dirac monopole is 
therefore: 

n 

Nevertheless, electrodynamics is perfectly consistent also without monopoles 
and therefore it does not require their existence. 

3.3.3 't Hooft-Polyakov monopoles 

It can be shown that a gauge theory coupled to scalar fields possesses monopole 
solutions, which are called 't Hooft-Polyakov monopoles [48]. 

The simplest example of non-abelian gauge theory having monopole solu- 
tions is the SO{3) Georgi and Glashow model [49], which is out by experiments. 
Nevertheless, the mathematical procedure developed for SO{3) can be applied 
to any compact non-abelian groups which break to a lower one containing some 
abelian factor, such as SU{b). 

Th SO{?>) lagrangian density is 



C = --FrF,^, + -D^^rD^r - V{<P) (S.SS) 



where 



{D^(j)r = d^r-ee'^^'Al^- 

V{ct>) = ^(^0.0-/^^ (3.89) 

This problem has a smooth solution with finite energy located in r = which 
behaves like a Dirac monopole at large distances since imposing total finite 
energy implies 

B~— (3.90) 

e r'^ 
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which has the same form of the last of fl3.78p . Therefore the field B is generated 
by a monopole having a magnetic charge 

9 = -- (3.91) 

e 

In this model the constant e is related to the electric charge operator by 

Q = eTa (3.92) 

where T3 is the third component of the weak isospin operatos; its smallest value 
is |. When e = | we get in modulus the Dirac condition fl3.87p with n = 1. 

It can be shown that stable monopole solutions occour for any gauge theories 
in which a simple group is broken down to a smaller group containing an explicit 
factor. This fact is compatible with the reasonable idea that electric charge 
quantization and the existence of monopoles are related and that the eletric 
charge quantization folows from the SSB of a simple gauge group. 



Experimental bounds 

Monopoles have many astrophysical and cosmological effects which depend on 
their average flux in the universe, defined as 

{Fm) = ^ (3.93) 

where hm is the monopoles number density and Pm stands for their velocity. 
No monopoles have been observed so far, therefore astrophysics and cosmology 
observations lead to only upper limits on (Fm)- 

The most stringent mass-independent limit on (Fm) is the Parker limit which 
is based upon the survival of the galactic magnetic fields: 

{Fm) < 10'^^ cm-^sr-hec^^ (3.94) 

This bound is analyzed considering also its dependence on the monopole mass 
in [50]. The larger possible flux of magnetic monopoles compatible with the 
survival of the galactic magnetic fields is 

{Fm) — 10"^^cm~^sr~^sec~^ 



for rriM ^ lO^^GeV, /?m ^ 3 ■ 10"^ (3.95) 

A direct experimental bound has been set by the MACRO collaboration [51]. 
The interaction of a magnetic monopole with a nucleon indeed can lead to a 
baryon-number violating process in which the nucleon decays into a lepton and 
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one or more mesons {catalysis of nucleon decay). Searching for these events in 
the MACRO detector lead to the bound: 

{Fm) < 3 ■ IQ-^^cm-^sr-^sec-^ 
for 1.1 ■ 10"^ < /9m < 5 ■ 10-^ (3.96) 



Chapter 4 

Milli-Charged Particles and the 
mirror universe 



The possible existence of particles having electric charge not multiple of g e has 
been longly investigated. These hypothetical particles as usually referred to as 
MCP, which stands for Milli- (or Mini-) Charged Particles. 

We have seen in chapter [3] that several advanced theories vinculate the elec- 
tric charges of the known elementary particles to be integer multiples of a fixed 
unity; nevetheless, it is possible to add new particles with unquantized charge 
to the model under some special conditions. 

For instance, Holdom [52] realized in 1985 that in a theory gauged by the 
group X U{1) there can be millicharged interactions; this scenario is com- 
patible with some gauge theories as well as with the constraints arising from 
the anomalies cancellation since the interaction takes place via a vector current 
which is automatically anomaly-free. 

We will explain in this chapter how the presence of a second group 
induces milli-charged interactions; then we will introduce a SM extension - called 
mirror universe - where this mechanism can be naturally introduced. 

4.1 Paraphotons and millicharges 

The mechanism by which millicharged particles are introduced is based on a 
lagrangian density containing two photons kinetically mixed. Two transforma- 
tions are performed, a non-unitary one to diagonalize the photon kinetic term 
and a unitary one to make ordinary matter coupled to one only photon. By mean 
of these transformations, particles which were coupled to the only non-ordinary 
photon get a small charge e with respect to the ordinary photon [52, 54] 

Let us start with a generic lagrangian density for a model gauged by the 
group W(l) xW(l) X .... It contains a mass term plus the interactions between 
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the photons and the fermions and a kinetic term which we assume to be non- 
diagonahzed: 

C = Ck^n + Cmass + Unt = -^F^MpF + ^A^ MaA +... (4.1) 

In the following we will only consider the photons' kinetic and mass term for 
sake of semplicity. As stated before, we will assume a diagonal mass matrix and 
a kinetic term having small quantities e out of the diagonal: 

M,= fl n °1 (4.2) 



V e 1 y ^ V m2 / 

With this choice of the matrices, the lagrangian in eq. (l4.ip contains the inter- 
action term for the photons: 

Ant = —F'^.F^'' (4.3) 

which is gauge invariant and renormalizable; it can exist at tree level and may 
be induced in GUTs [56,64]. 

To get a canonical kinetic term, we transform the kinetic matrix Mp to a 
diagonal matrix D via: 

Mp = ODO^ 0^0 = 1 (4.4) 

which leads to a canonical kinetic term under the non-unitary transformation 
of the fields: 

F' = D^20^F A' = d'W^A (4.5) 

Applying these transformations the lagrangian density becomes: 

C = -^F'^h^.F' + \a'^M'^A' (4.6) 



where 



1 



M'^ = D—^O^MaOD—^ (4.7) 
An explicit working out of these matrices gives: 

Let us now insert the physics in the model. Since the ordinary photon is 
massless, we assume that both Ai and A2 are massless and therefore state 

mi = 7712 = 0. 

Any pair of orthogonal combinations of the A'^ fields will keep canonical 
kinetic term. The theory can contain particles coupled to one only photon or 
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both, but since QED has been tested to very high precision and its previsions 
are in very good agreement with the experimental results, we want ordinary 
particles - electron, muon and so on - to couple to one only photon. We therefore 
re-define the fields such that one photon - which we will call Ai - will satisfy 
this condition. The field Ai therefore represents ordinary QED photon, while 
the orthonormal combination A2 will be called paraphoton and will be sterile 
with respect to the ordinary matter. 

If the model contains ordinary particles coupled to A^ with charge e[ - which 
we will call /i - and para-particles coupled to A2 with charge e'g - which we will 
call /2 - and the interaction lagrangian density is of the usual form: 

Cint = -e;/i7^/iA5' - e'j2ij2A^ (4.9) 

performing the transformations we introduced above we get: 

• ordinary particles /i couple to the only Ai with the usual electric charge 
Ci oc e[] as a consequence the electric charges of ordinary particles keep 
the same ratio they have in the SM scenario and - for istance - the electron 
to proton charge ratio is not modified; 

• para-particles /2 couple to A2 with charges 62 oc e'2 but also to the ordinary 
photon Ai with charge e. 

In conclusion: paraparticles, which were coupled to the only paraphoton in the 
lagrangian (14. ip . via this mechanism gain a small electric charge which makes 
them coupled also to ordinary photon. This small charge is equal to e, which is 
the kinetic mixing parameter. 

4.2 The mirror universe 

We introduced in §1.41 that elementary particles experience left-handed weak 
interactions. But it has been suggested since 1956 by Lee and Yang [17] that 
there may exist right-handed particles restoring the left-right symmetry. Ten 
years later Okun and Pomeranchuk [18] suggested that these particle may come 
from a hidden mirror sector, which is an exact copy of our one but where weak 
interactios are right-handed. 

The mirror theory is based on the gauge group G x G' , where ' marks quan- 
tities of the mirror sector. If it is not differently specified, we will assume below 
that the discrete symmetry G ^ G' , which interchanges corresponding fields of 
G and G' and is called mirror parity, is valid; this symmetry guarantees that 
two particle sectors have identical lagrangians with the same interactions and 
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coupling constants. If the electro-weak symmetry scale is different in the two 
worlds, this parity will be spontaneously broken and this can lead to different 
physics in the mirror sector [66,71-73]. 

The two worlds communicate through the gravity, but there may also be 
mini-charged interactions like the ones we introduced in §4.11 Other possible 
interactions between the two worlds have been investigated, for instance the 
mixing of ordinary and mirror (sterile) neutrinos [66,74] or ordinary and mirror 
neutron [63], a quartic interaction between the two worlds' Higgses [62,73] and 
so on - see [62] for more references. 

An interesting feature of the millicharged mirror world is that every particle 
contained in the SM has its mirror counterpart, which has, with respect to 
the ordinary photon, the same electric charge of its standard corresponding 
but scaled by a factor e. Hence, the milli-charged mirror sector can not give 
rise to exotic unwelcome phenomena, such as violation of the electric charge 
conservation or atoms having the same number of mirror protons and electrons 
but nevertheless globally charged. Moreover, the conditions which ensure the 
anomaly cancellation in the SM - see §3.1.51 - if satisfied in the ordinary sector, 
will be automatically satisfied in the whole theory. 

Mirror standard model and mirror symmetry 

Let us consider the mirror principle applied to the SM; the whole gauge group 
for this case will be G x G" = SU{3) x SU{2) x U{1) x 5W(3)' x 5W(2)' x 
Also the lagrangian is the same than in the SM; the mirror symmetry 
introduced above, which interchanges all ordinary particles with their mirror 
partners, induces also the lagrangian interchange 

C^C (4.10) 

such that the total lagrangian C + C remains invariant. 

The SM particles content has been already introduced in §1.51 - see equations 
f ll.45p . fll.46p and fll.47p . The mirror sector contains the same particles, which 
will be marked by the usual '; also the barion and lepton numbers B',L', L[ i = 
e, /i, T are assigned in analogy with the ordinary ones. 

When we introduced the SM in chapter [1] we stressed that weak interactions 
are not symmetric between the left and the right components: parity P : ^ 
/r and CP : fi ^ ffj simmetries are both broken 

Nevertheless, it is possible to preserve the symmetry between left- and right- 
handed components in the mirror scenario. This symmetry - sometimes called 
matter parity - is made through the transformations [53] : 
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/c 



^ = 1,2,3 , j = l,...,8 (4.11) 

It is worth to stress that if we require the vahdity of this symmetry, the intro- 
duction of the mirror world does not add any new parameter to the SM; the 
matter symmetry vahdity indeed also imphes that the gauge couplings and the 
Higgs potential are the same for the two sectors, while the Yukawa couplings 
must satisfy the property: 

n,.,. = n',V (4.12) 

Through these transformations, the left-handed particle physics we have in the 
standard model has an identical but right-handed counterpart in the mirror 
sector. The left-right symmetry is therefore restored in the two sectors theory. 

Generalizations to supersymmetric and grand unified models can be found 
in literature - see for instance [62,64,65]. From grand unified mirror models 
we can get a natural understanding of the smallness of the mixing parameter e. 
Indeed, cosmology, astrophysics and experiments bound e to have small values 
(of order < 10~^, see §4.4p : but if the theory does not provide any reasons for 
this smallness, it would be natural to expect e ~ 1. 

A kinetic mixing term like (14. 3 p is forbidden in GUTs, like for instance 
SU{5) X SU{5y which do not contain abelian factors [65]. However, given that 
both SU{5) and SU (5)' symmetries are broken down to their SU (3) x SU (2) x 
U{1) subgroups by the Higgs 24-plets $ and <!>', such a term could emerge from 
the higher order effective operator 

C = -^iG^-<t>)iG'^,<t>') (4.13) 

where Gfj^^ and G'^^ are field-strength tensors respectively of SU{5) and SU{5)', 
and M is some cutoff scale which can be of the order of Mp or so. The operator 
(Km can be effectively induced with e ~ C{{^)/Mf {( ~ a/Sn ~ lO'^) by 
loop-effects involving some heavy fermion or scalar fields in the mixed repre- 
sentations of SU{5) X SU{5)'. It can be shown that, taking the GUT scale as 
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($) 1016 GeV and M ~ Mp, the strength of kinetic mixing term (14.31) . could 
vary vary from e ~ 10^^° to 10^^. 

Finally there is a last remark: also if the micro-physics is the same but 
with opposite chiralities in the two sectors, these sectors may have different 
macroscopic evolutions. In particular, as we will see more in detail in §4.3[ 
BBN bounds the mirror sector to a temperature lower than ordinary world. 



4.3 Mirror cosmology: thermodynamics 

We introduced in §2.51 the criterion through which, by comparing experimen- 
tal measures and theoretical calculations of the primordial abundance of light 
elements, the number of extra neutrinos is bounded to 

AN^ < 1 at 95% C.L. (4.14) 

Since the mirror universe particle content would lead to 

AAT. = ^ = 1^ ^ 6.14 (4.15) 
2- I 1.75 ^ ' 

it is straightforward that the mirror sector can exist only if there is some mech- 
anism which reduces this number. This is naturally achieved if the mirror world 
temperature is lower than the ordinary one - see equation (12.411) : in particular 
equation (I4.14p is satisfied if 

T' 

— < 0.64 (4.16) 

As stated above, it is possible that the two sectors, in spite of having the same 
microphysics, may have different macroscopic evolutions. In effect, the mirror 
sector temperature can be lower than the ordinary one if: 

• The mirror sector has the inflationary reheating temperature lower than 
the ordinary one, as it happens in certain models [68-70]; 

• Interactions between ordinary and mirror particles are so "weak" that the 
worlds do not come in thermal equilibrium. 

The expressions for entropy and energy density we introduced in §2.3.11 can be 
easily generalized to the mirror world [84]: 

s'iT') = '^-^qUT)r' (4.17) 
qUT') ^ E ^'biT')(§X + 1 E ^'fiT')(^X (4.18) 

bosons \ / fermions \ / 
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p'(T') 



TC 

30 



/4 



(4.19) 



bosons \ / fermions \ / 

If the two sectors communicate only through gravity and they both expand 
adiabatically, they evolve with separately conserved entropies and therefore the 
ratio X, defined as: 



X = 



T'{t) 



X 



q'{T') 



(4.21) 



V s J m 

is time independently]. 

The universe evolution in the Friedmann equations - see fl2.10p - is deter- 
mined by the total energy density p = p + p' ., which enters the Hubble constant 
as 



H 



which can be recasted using the definitions fl2.40p and (14.191) as 



H = lMJg{T) 



Mf 



l.mJ-g'{T'\ 



■T 



12 



where 



giT) = g{T){l + ax') g'{T') = g'{T') 1 + — 

V ax^ 

and the factor a{T,T'), defined in the following way 



(4.22) 



(4.23) 



(4.24) 



9'iT') 

giT) 



q'{T') 



(4.25) 



takes into account that for T ^ T' the relativistic particle content of the two 
worlds can be different. 

It can be easily seen that the bound (14.141) on the extra-neutrino number 
can be expressed as a bound on x. If we assume indeed ~ in (14.211) and 
recall the condition on the mirror temperature in (14.161) we have 



X 



r 



< 0.64 



(4.26) 



"'^ Other interactions may lead to entropy exchange between the two sectors. This is analyzed 
more in detail in chapter [6] for minicharged interactions. Of course entropy exchange imply 
that the parameter x is not a constant anymore but a function of time (or temperatures). 
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The link between x and AiVj, can be summerized in the formula 

AN^ = 6.Ux^ (4.27) 

which is valid in the limit a ^ 1. This means that the bound on x we reported in 
equation fl4.26p has a very weak dependence on AA^^, but also that any changes 
on X has a great impact on AN^. Some consequences of this phenomenon are 
analyzed in §6.2[ 

In conclusion, the presence of the mirror sector having x < 0.64 has pratically 
no impact on the standard BBN. But the impact of the ordinary world for the 
mirror BBN is much stronger, as we will discuss further in the following chapter. 

4.4 Mirror matter as dark matter 

Cosmological observations indicate that the Universe is nearly fiat, with the 
energy density very close to the critical one: Qtot = 1- The non-relativistic 
matter in the Universe consists of two components, baryonic (B) and dark (D). 
Dark Matter (DM) is non-luminous and non-absorbing matter which interacts 
with the ordinary one mainly via gravity; its existence is now well established, 
since it has measurable effects on the velocity of various luminous objects and 
on the cosmic microwave background anisotropics [86]. 
Recent data fits imply [38]: 

^3/^2 = 0.0222 ± 0.0007 f^o/^^ = 0.105 ± 0.004 (4.28) 

where h = 0.73 ± 0.02 is the Hubble parameter. Hence, matter gives only a 
fraction of the total energy density: 

nM = nB + nB = 0.24 ± 0.02 (4.29) 

while the rest is attributed to dark energy (cosmological term): 

= 0.76 T 0.02 (4.30) 

The closeness of and (^d/^b = 4.7 ±0.3) gives rise to a painful problem, 
called Fine Tuning problem. Both pD and pe scale as ~ with the Universe 
expansion, and their ratio should not be dependent on time. Why then these 
two fractions are comparable, if they have a drastically different nature and 
different origin? 

An answer to this question is found assuming that DM is made of mirror 
baryons. We can indeed define the baryon mass density p-Q = m^n^, where 
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TUB — 1 GeV is the nucleoli mass, and ub is the baryon number density. Thus 
we have 

l^B/i' = 2.6 X 10^ — (4.31) 

s 

where s is as usual the entropy density; equation (14.281) translates in B = ub/s ^ 
0.85 X 10"^°, in a nice consistence with the BBN bound B = (0.5-1) x 10"^° [86]. 

The origin of non-zero baryon asymmetry B, which presumably was pro- 
duced in a very early universe as a tiny difference ub = — ni between the 
baryon and anti-baryon abundances, is yet unclear. The popular mechanisms 
known as GUT Baryogenesis, Leptogenesis, Electroweak Bariogenesis, etc., all 
are conceptually based on out-of-equilibrium processes violating B{B — L) and 
C/CP and they generically predict i? as a function of the relevant interaction 
strengths and CP-vioUating phases. 

Concerning dark matter, almost nothing is known besides the fact it must be 
constituted by some cold relics with mass mo which exhibits enormous spread 
between different popular candidates as e.g. axion (~ 10~^ eV), sterile neutrino 
(~ 10 keV), WIMP/LSP (~ 1 TeV), or Wimpzilla (~ 10^^ GeV). None of these 
candidates has any organic link with any of the popular baryogenesis schemes. 
The respective abundances rij) could be produced thermally (e.g. freezing-out of 
WIMPs) or non-thermally (e.g. axion condensation or gravitational preheating 
for Wimpzillas), but in no case they are related to the CP- violating physics. In 
this view, the conspiracy between = fnunn and pb = itibUb indeed looks as 
a big paradox. 

Mirror baryons may shed a new light to the baryon and dark matter coinci- 
dence problem: they are invisible in terms of the ordinary photons but interact 
gravitationally with ordinary matter and thus they constitute a viable dark 
matter candidate [68,71,84]. 

In general we can write 

= Q!^ + n'^ < (4.32) 

where Q'^ is the mirror baryon density and is some kind of exotic dark 
matter. We are interested in a situation when the ratio 

f3 = ^ (4.33) 

falls in the range from 1 to few tens. From (3 we can also work out a lower 
bound on the parameter x: since n'^ = x^n^, we can write 

^ = x^"^ ^ x^=(3^ (4.34) 
r] rj' 
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Recalling now that ~ 10 ^ and 77' < 10 ^, we get the lower limit 

X > 10"^ (4.35) 

Let us assume now Qb = ^d] the closeness of 17/) = Qb and Q'^ arises because of 
the mirror symmetry. Their difference instead can be explained, in a symmetric 
mirror scenario where tub = m'^, in terms of a n'^ > n^, as one can achieve in 
certain leptogenesis models - see [62,64,85]. 

An other intriguing mechanism leading to Q'b > I^b is based on the possi- 
bility that the Higgs Vacuum Expectation Values (VEVs), which is the electro- 
weak scale of the theory, may be different in the two sectors [68] : 

{<P)=v {<f>')=v' ^ = C^1 (4.36) 

Changing the Higgs VEV implies different weak interaction strenghts and dif- 
ferent masses of the weak gauge bosons and of the fermions in the two sectors: 

G'f = ^ M{^^z = (Mw,z mf = Cm'j (4.37) 

and also a slower change in the confinement scale A, which scales approximately 
as [71,72]. This way, assuming ( ~ 100 we find A' ~ 800Me^, which implies 

m'j ~ 100 m'^ 5mB (4.38) 

This discrepancy on the baryon mass leads to Q'b ^ 5Qb [67]. Consequences 
on BBN of having heavier mirror fermions will be analyzed in section [6l Other 
effects arise since the Bohr radius scales as ~ 0.01 and this implies that 
mirror atoms are much more compact and therefore less collisional than the 
ordinary ones. Moreover, the hydrogen recombination and photon decoupling 
in M-sector would occur much earlier the matter-radiation equality epoch, and 
as a consequence, mirror matter will manifest rather like a cold dark matter. 

Also, mirror nuclei will be all, but hydrogen, unstable since the light quark 
mass difference scales as 

(m^ - m'J ^ ({ma - m„) (4.39) 

We therefore expect the mirror neutron n' to be heavier than the mirror proton 
p' by few hundred MeV. Clearly, such a large mass difference cannot be com- 
pensated by the nuclear binding energy and hence even bound neutrons will be 
unstable against (3 decay n' p'e'u'^. Thus in such an asymmetric mirror world 
hydrogen will be the only stable nucleus [68]. 

In the following, if not differently specified, we will always refer to the sym- 
metric mirror model, that is we will assume equal masses for the elementary 
particles. In this scenario mirror matter can be unambigously described by 
mean of two only free parameters, x and f3. 
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4.5 Bounds on MCPs and prospects for their 
search 

Bounds on MCP parameters - mainly charge and mass - have been worked out 
by several authors - see for instance [91]. Complete and more recent reviews can 
be found in [92] and [94], where bounds arising from accelerator experiments 
(AC), BBN, globular clusters (GC), supernova 1987A (SN) and white dwarfs 
(WD) are reported. Also, Dubovsky et al. derived bounds from the Cosmic 
Microwave Background Radiation (CMBR) [93]. 

We are interested in a scenario where the millicharged particles have mass 
~ TUe = 5 ■ lO^eV^ or ~ 100 rrie = 5 ■ 10''eV^. The e values ruled out for 
these masses in the previous references are listed in Table 14. 1[ The BBN bound 



Table 4.1: BBN bounds fom literature. 





BBN 


SN 


AC 


~ rrif. 


e < 2 ■ 10-9 


10-^ < e < 10-9 


e < 10-5 


~ 100 me 






e < 10-=^ 



is for Davidson et al. [92]: 

e < 2 • 10-9 vahd when m, < rUe (4.40) 

This bound has been calculated for a model containing only MCPs and one 
paraphoton in the ultrarelativistic (UR) limit, that is assuming massless MCPs 
and electrons and approximating the interaction rate F to - see §6.1. It 

ruR = |^(ea)'T (4.41) 

But BBN takes place at temperatures comparable with the electron mass - 
see §2.51 - when this approximation is not very goocj^ as we show in Figure 14. 1^ 
where the ratio 

r ratio = (4.42) 

TuR 

is plotted. 

The bounds we listed above have been calculated for models containing 
MCPs with or without one paraphoton but are not specific for the millicharged 
mirror model. Glashow pointed out [56] that a sensitive laboratory test of 



^The procedure we used to numerically work out T can be found in §6.1.11 
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Figure 4.1: Ratio ofT numerically calculated and T in ultrarelativis- 
tic particles approximation. 

the mirror model comes from the orthopositronium system: the kinetic mixing 
leads orthopositronium (o-Ps) to oscillate into mirror orthopositronium with 
the probability: 

P{o - Ps ^ o - Ps') = sin^ Ljt, (4.43) 

where u = 27re/ and where / = 8.7 x 10^ MHz is the contribution to the 
ortho-para splitting from the one photon annihilation diagram involving or- 
thopositronium. Since mirror o-Ps are not detected, in experiments they look 
like invisible decays. Recent searches by Badertscher et al. gave no evidence for 
positronium decays into invisible states in a collected sample of about 6 ■ 10^ 
decays; this way the photon mixing strenght is bounded to e < 1.55 ■ 10^'' at 
90% C.L. [57]. 

BBN bounds specific for the mirror model have been worked out in 1987 
by Carlson and Glashow [55], who derived e < 3 ■ 10^®. Nevetheless, we have 
today new data and new theoretical hints which render this limit worth to be 
reanalyzed. This will be done in chapter [61 

An other interesting experiment which may have detected millicharged in- 
teractions of ordinary and mirror particles is DAMA/Nal [60,61]. DAMA is an 
observatory for rare processes based on the development and use of various kinds 
of radiopure scintillators located in the underground laboratory of Gran Sasso. 
The DAMA/Nal experiment was searching for scattering of WIMPs (weakly 
interacting massive particles, which are a popular DM candidate) on a target 
of 100 kg of radiopure Nal, measuring the recoil energy of the scattered Nal 
atoms. The detected signal is expected to have a small annual modulation due 
to the Earth's motion around the sun: the Earth should be crossed by a larger 
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flux of Dark Matter particles roughly around June 2nd (when its rotational ve- 
locity is summed to the one of the solar system with respect to the Galaxy) 
and by a smaller one roughly around December 2nd (when the two velocities 
are subtracted), according to 



The experiment took data over 7 annual cycles and was put out of operation in 
July 2002. The data analysis gave evidence for an annual modulation at 6.3a 
CL. with the fitted parameter values T = (1.00 ±0.01) year (the expected value 
for T was 1 year), to = 144 ± 22 days (the expected value for to was 152 days, 
that is 2 June), while the signal strength is A = (0.019 ± 0.003) cpd/kg/keV. 

Foot proposed to interpret these data in terms of elastic scattering of mil- 
licharged mirror atoms on ordinary matter nuclei in the symmetric mirror sce- 
nario [58,59]. Indeed, he showed that DAMA/Nal as well as other DM searching 
experiments (such as CRESST/Sapphire) are sensitive to this kind of interac- 
tions provided that e ~ 5 ■ 10~^. We will show in chapter [H] that such a value 
for e in a symmetric mirror model, where the particle masses are the same than 
the ordinary ones, is compatible with BBN bounds, also if it may give prob- 
lems when considering CMB and LSS. We will then state reasons for further 
investigations of the same phenomenon but within an asymmetric mirror model, 
having m'j ~ 100m/, m'^ ~ mjfl. 




(4.44) 



■^Where / stands for any fermions and B for any baryons. 



Chapter 5 

BBN in presence of the mirror 
world 

BBN in presence of the symmetric mirror sector, where the particle masses are 
the same than the ordinary ones, is analyzed in this section. In §5.1l we calculate 
the equations which link the ordinary and the mirror sector temperatures; then 
we numerically solve these equations using Mathematica. Once the tempera- 
tures are known, it is possible to work out the total exact number of degrees of 
freedom (DOFs) in both sectors, which can be, as it is common in the literature, 
expressed in terms of extra-neutrino number. We report in §5.1.21 the results 
of these calculations, which have been made again with Mathematica. Finally, 
we use the Kawano code for BBN [87] to work out the mass fraction of light 
elements produced in both sectors during the primordial nucleosynthesis; the 
presence of the other sector indeed, leads in both sectors to the same effects 
of having more particles. These topics can be found in §5.31 and §5.4[ For the 
moment we do not take in account of the possible existence of MCPs or other 
interactions - but gravity - between the two secors; millicharged interactions 
will be instead considered in the following chapter. 

5.1 Neutrino decoupling and annihilation: 
working out the mirror and ordinary tem- 
peratures 

As already stressed in chapter HI ordinary and mirror sectors have the same 
microphysics; it is therefore obvious that the neutrino decoupling temperature 
Tdj, is the same in both them, that is T^i, = T'^^. 

This fact will be used together with the entropy density conservation to 
find equations the solution of which will give the mirror photon temperature T' 
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corresponding to any values of the ordinary photon one T. 
5.1.1 The equations 

We argumented in §4.31 that the mirror world must be colder than the ordinary 
one and therefore the neutrino decoupling takes place before in the mirror sector. 
If we call Td;/' the ordinary world temperature when the mirror annihilation 
takes place, we can split the universe evolution in three phases: 

Phase 1: T > Td^, 

Photons and neutrinos are in thermal equilibrium in both worlds, that is T^, = 
T , Tl = T'. Using equations (l238|) . (IXBOj) . (ICTj) and fICTD and the corre- 
sponding ones for the mirror sector that is fl4.17p . fl4.18p . (14. 19p and (I4.20p we 
are able to calculate the DOFs number in ordinary or mirror world; but to work 
out the total DOFs number, that is summed on both worlds, we need the mirror 
temperature T' as a function of the ordinary temperature T or viceversa. 

If we neglet the entropy exchanges between the sectors we can get both these 
functions imposing x =const in: 



o s ■ 

x-^ = 

s ■ a-^ 



he{T') + q^ + hv 



rp/3 

{5.1] 



T3 



_8?e(T) + qj + ^q, 

where qu = G and = 2 as introduced in (I2.42p . while qe(T) stands for 



T 

(5.2) 



45 



and Pe{T), PeiT) are the same than in (I2.34p and (I2.35p . 

Equation (15. ip can be numerically resolved and will finally give the function 
T'(T) for every T > T^u'- 



Phase 2: Tdu<T < Tdu' 

At T ~ Td^' mirror neutrinos decouple - see §2.41 - and soon after mirror elec- 
trons and positrons annihilate, raising the mirror photon temperature. Never- 
theless, ordinary photons and neutrinos still have the same temperature T. We 
have therefore two equations: the first one is 
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and its solution gives as a function of T'; the second is 



X 



s ■ a-" 



|ge(T') + 


rplZ 1 7 

8 


T' 3 




l<le{T) + q^ + 


y3 



(5.4) 



and its solution gives T' function of T. 
Phase 3: T < T^^ 

At T ~ Tciy ordinary neutrinos decouple and soon after ordinary electrons 
and positrons anihilate; therefore we need one more equation to work out the 
ordinary neutrino temperature as a function of the ordinary photon one T: 



22 
21 

22 
21 



qe{T') + q, fr^ 
' \ T' 



heiT) + q, (T 



T 



(5.5) 



(5.6) 



X 



s' ■ a? 



s ■ a-" 



heiT') + q^ 


T" + lq.K^ 


lqe{T) + q^ 





(5.7) 



Once both ordinary and mirror photon temperatures are known, it is straight- 
forward to calculate the total energy and entropy densities; then, reversing 
equations (12.381) and fl2.40p . we can work out the entropic (g) and energetic {g) 
total number of degrees of freedom. Calculations have been made for several 
different values of reported in the following section. 

5.1.2 Numerical solution 

The equations we introduced above have been numerically solved with Mathe- 
matica. Since x is a free parameter in our theory, several values have been used 
for it - from 0.1 to 0.7 with step 0.1 or less, see §4.31 - to study its effect on the 
light elements production. In the extreme cases 



T > Tdu' 



T, 



Du 



X 



or 



IMeV 



(5.8) 



we expect q{T) ~ q'(T'); therefore in this limit - see fl4.2ip - the ratio of mirror 

len 

(5.9) 



and ordinary photon temperature should be x, that is ^ = 1. Instead, when 
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x=0 . 7 

x=0 . 6 

x=0 . 5 

x=0 . 4 

x=0 . 3 

x=0.2 

x=0 . 1 



T 

Figure 5.1: ^ for several values of x. The asymptotic value for this 
ratio is 1 for both high and low temperatures, as expected. 

we expect q'{T') < q{T) because the annihilation takes place before in the 
mirror world, leading to a decrease of q'{T') and a corresponding increase of T' 
in order to keep constant the entropy densities ratio; thus we expect ^ > 1. 
Later on, when the ordinary electrons-positrons annihilate, they make even T 
increase and thus the ratio ^ decrease to the asympotic value 1. These remarks 
have been numerically verified; the ratio ^ is plotted in Figure 15.11 

5.2 The number of degrees of freedom 

Using equations from fl2.34p to fl2.36p to reverse (12.381) and (12. 40 p . we can work 
out the total number of entropic (g) and energetic {g) DOFs at any temperatures 
T. Of course, we can apply the same procedure to work out the standard, that 
is in absence of the mirror world, total q and g, as well as q' and g' and the 
influence of a world with the other one. 

In general we expect q (g) to have a cubic (quartic) dependence on x - see 
equation (I4.24p - when the temperature is not close to the u decoupling and the 

annihilation phases, that is, when a = 1 - see equation (14.251) . 

Numerical calculations have been made using Mathematica; in Table 15.11 
some values are reported for special temperatures and several values of x. As 
expected, the DOFs number is always higher than the standard and increases 
with X. Moreover, the mirror sector values are higher than the ordinary ones of 
a factor of order x~^ (for q) or x"^ (for g). 



T mirror 

X T Tmirror/xT 

1 . 5 I ' ' ' ' 1 ' ' ' ' — I ' ' ' ' — 



1.4- 



1.3- 



1.2- 



1.1- 




12 3 4 
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In Figure 15.21 q, q, q' and the ordinary DOFs number increase due to the 
mirror sector are plotted for the intermediate value x = 0.5; the second plot 
in the same figure shows the corresponding gs. In the figure q' {g') has been 
scaled by a factor {x^)', this way the asymptotic values are the same than the 
ordinary sector ones because at the extremes ^ = x (see Figure [5TT1) . We can 
see that q, g, q' and g' have similar trends, but, due to T' < T, mirror DOFs 
begin to decrease before. 

In Figure [HT5] q in ordinary and mirror sectors are plotted in comparison with 
the standard for several x values (from 0.1 to 0.7 with step 0.1). The foreseen 
quartic dependence of 5^ on x at the extremes is proved correct. We note that 
the plot shape does not change with x in the ordinary sector, while it does in the 
mirror one. This sector evolves with the temperature T' ~ xT < T; therefore 
the number of DOFs converges later for lower xs. 

The number of neutrinos in the ordinary sector 

We know from chapter [T] that the SM contains three neutrino species; the possi- 
ble existence of a fourth neutrino has been longly investigated, also using BBN 
constraints. This is why in the literature one can often find bounds on the num- 
ber of DOFs in terms of extra-neutrino number AN,^ - see for instance [77,89]. 
In general, the number of neutrinos N^, is found assuming that all particles con- 
tributing to the universe energy density, to the exclusion of electrons, positrons 
and photons, are neutrinos; in formula that means: 

g{T) = g^^(T) + g, + 'l-2K-(^ 

We can also work out the number of extra-neutrinos AA'^i,, that is the number 
of non standard DOFs in term of addictional neutrinos: 

g -Jhtd fT_" ^ 

■ 2 ■ It. 



AiV. = V^F-l^^r) (5.11) 



Both and AN^ have been worked out with Mathematica; some data are 
reported in Table 15.21 while plots for several x values and any temperatures 
from to 3 MeV are shown in Figure 15. 4[ We stress that the standard values 
Njj = 3 and AN^, = are the same at any temperatures, while a distinctive 
feature of the mirror scenario is that the number of neutrinos raises with the 
temperature. Anyway, this effect is not a problem; on the contrary it may be 
useful since some recent data fits give indications for a number of neutrinos at 
recent times higher than at BBN [89]. 
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Table 5.1: Standard and not standard total DOFs number for several 
X values in both ordinary and mirror sector. 



T{MeV) 


Ti = 0.005 


T2 = 0.1 


T3 = 0.5 


T4 = l 


T5 = 5 


q (standard) 


3.91 


4.78 


10.0 


10.6 


10.74 


g (standard) 


3.36 


4.30 


10.0 


10.6 


10.74 


X = 0.1 


T'{MeV) 


0.0005 


0.0107 


0.0675 


0.124 


0.511 


q 


3.91 


4.78 


10.0 


10.6 


10.75 


q' 


3910 


3910 


4070 


5520 


10000 


9 


3.36 


4.30 


10.0 


10.6 


10.74 


9' 


33600 


32900 


30000 


44400 


98400 


x = 0.3 


T'{MeV) 


0.0015 


0.0321 


0.170 


0.315 


1.50 


q 


4.015 


4.91 


10.3 


10.8 


11.0 


q' 


148 


148 


261 


347 


404 


9 


3.39 


4.34 


10.1 


10.6 


10.8 


9' 


418 


409 


750 


1080 


1320 


X = 0.5 


T'{MeV) 


0.0025 


0.0533 


0.263 


0.508 


2.50 


q 


4.40 


5.38 


11.3 


11.9 


12.1 


q' 


35.2 


35.5 


77.3 


90.5 


96.5 


9 


3.57 


4.57 


10.7 


11.2 


11.4 


9' 


57.2 


56.6 


138 


168 


182 


x = 0.7 


T\MeV) 


0.0035 


0.0733 


0.357 


0.704 


3.50 


q 


5.25 


6.42 


13.5 


14.2 


14.4 


q' 


15.3 


16.3 


4.88 


40.6 


42.0 


9 


4.17 


5.35 


12.4 


13.1 


13.3 


9' 


17.4 


18.5 


47.8 


53.3 


55.5 
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Energy DOFs (x=0.5> 




Figure 5.2: In this figure we plot: standard q and g (in dark green), 
mirror contribution to total ordinary q and g (in light green), total 
q and g in mirror (in blue) and ordinary (in red) sectors, all for 
X = 0.5. Mirror q andg (blue lines) have been multiplied respectively 
by and x^ to make them comparable with the ordinary values since 

T' 1 

^ ~ i. 




Figure 5.3: g in ordinary and mirror sectors for several values of x 
compared with the standard. 
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Table 5.2: The number of neutrinos in the ordinary sector for some 
special cases. 



T{MeV) 


standard 


X = 0.1 


X = 0.3 


X = 0.5 


X = 0.7 


ordinary sector 


0.005 


3.00000 


3.00074 


3.05997 


3.46270 


4.77751 


0.5 


3.00000 


3.00074 


3.05997 


3.40706 


4.52942 


1 


3.00000 


3. .00071 


3.05202 


3.39166 


4.49133 


5 


3.00000 


3.00063 


3.04989 


3.38430 


4.47563 



TOTAL NEUTRINO NUMBER 



ANEUTRINO NUMBER 



1.75 
1 . 5 
1 . 25 

1 

0.75 
. 5 



x = 


7 


x = 


6 


x = 


5 


x = 


1 


x = 


3 


x = 


2 


x = 


1 



Figure 5.4: N^, and ANi, in the ordinary sector for several x values. 
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Predictions for special temperatures It is possible to work out and 
AA^,^ at T and T -C Tanne^ starting from the standard values - see (12.431) 

and flCTjl : 

(7(r > T^,) = 10.75 (?(r<r,„„e±) ^3.36 (5.12) 

Without the mirror universe, we have, as expected: 



10.75-2-1-4 

N,{T:^Td..) = ' =3 



3.36-2 /ll\l 



N,{T<^Tanne±) = \ . " ( — ) =3 (5-13) 



• 2 V 4 / 



While, when the mirror universe is present, we have g = g{l + x^) at the special 
temperatures we are considering; hence 



10.75(1 + -2- I -4 



^ X 3.36(1 + 2 /11\3 ^ 

N,iT<^Tanne±) = ^7^^ ) (5.14) 



8 

From the equations written above we can see that the rise in is: 

■1.25-x^ (5.15) 



4 1 

X 



■ 2 



10.75 - 3.36 ^ ^ 



11\ 3 
T 



This leads to Ni, higher than the standard in presence of the mirror world and 
to a further rise of this parameter at low temperatures. For instance, if x = 0.7 
equations (15.141) give: 

iV,(T > T^,) ~ 4.5 and Ar,(T < T,„„e±) ^ 4.8 (5.16) 
The number of neutrinos in the mirror sector 

Similarly the number of neutrinos in the mirror sector can be worked out as: 



1-2 [t: 



Once again these values are higher than the ordinary ones by a factor x^; they 
have been numerically computed using Mathematica and some special values 
are reported in Table 15.31 
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Table 5.3: The number of neutrinos in the mirror sector for some 
special cases. 



r (MeV) 


X = 0.1 


X = 0.3 


X = 0.5 


X = 0.7 


mirror sector 


0.005 


74011 


917.0 


121.4 


33.83 


0.1 


62007 


805.0 


111.6 


32.21 


0.5 


61447 


763.1 


101.9 


28.86 


1 


61435 


761.8 


101.4 


28.66 


5 


61432 


761.4 


101.3 


28.59 



5.3 The light elements production in the ordi- 
nary sector 

As we have seen in §5.21 the presence of the mirror sector can be parametrized 
in terms of extra DOFs number or extra neutrino families; therefore, since the 
physical processes involved in BBN are not affected by the mirror sector, we 
can use and modify a pre-existing BBN code to work out the light elements 
production. 

We chosed the Kawano code for BBN [87] since it is a well-tested and fast 
program and its accuracy is large enough for our purposes. The Kawano code 
is an updated and modified version of the code of R.V. Wagoner [88]. It solves 
for elemental abundances arising from the epoch of primordial nucleosynthe- 
sis in the early universe by mean of a second-order Runge-Kutta driver. The 
parameter values we used are 

ri = 6.140 • 10~^° T = 885.7 s (5.18) 

where rj is the final baryon to photon ratio. The number of DOFs enters the 
program in terms of neutrino species number; this quantity is a free parameter 
and can set by the user. But the program uses the same number during the 
whole BBN process. 

Instead, the mirror sector leads to variable A^;^! We therefore used before 
a first order approximation where AN^, = const oc (the so called zero ap- 
proximation) to get a feehng of the nuclides production. Then, in order to get 
more accurate results, we modified the program to read the neutrino number 
for any temperatures and x values as an external input. This input comes from 
a file, which has been written with a Mathematica program which solves the 
equations which describe the thermodynamical evolution of both sectors - see 
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Figure 5.5: The development of primordial nucleosynthesis for the 
standard universe (that is without the mirror sector) and in the two 
sectors when the mirror is present. To work out these data we used 
the parameters a; = 0.7, /5 = 5 . 
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§5.11 and §5.2[ The only parameter of the mirror sector which affects ordinary 
BBN is x; the baryon ratio P - defined in equation fl4.33p - does not induce any 
changes on the production of ordinary nuchdes. We will see in the following 
section that (3 plays a crucial role for the mirror nuclides production. 



Table 5.4: Light elements produced in the ordinary sector. 



elements 


std 


X = 0.1 


x = 0.3 


X = 0.5 


X = 0.7 


n/H {lO-'"^) 


1.161 


1.161 


1.159 


1.505 


2.044 


P 


0.7518 


0.7518 


0.7511 


0.7463 


0.7326 


D/H (10-^) 


2.554 


2.555 


2.575 


2.709 


3.144 


T/H (10-^) 


8.064 


8.065 


8.132 


8.588 


10.07 


^He/H (10-5) 


1.038 


1.038 


1.041 


1.058 


1.113 




0.2483 


0.2483 


0.2491 


0.2538 


0.2675 


(10-14) 


1.111 


1.111 


1.124 


1.210 


1.499 


'^Li/H (10-^^) 


4.549 


4.548 


4.523 


4.356 


3.871 


^Be/H (10-10) 


4.266 


4.266 


4.238 


4.051 


3.502 


«Li + /H (10-15) 


1.242 


1.242 


1.251 


1.306 


1.464 



In Table |531 the program output, that is the light elements produced in the 
ordinary sector at the end of BBN process (at T ~ 8 ■ 10"^ MeV), is reported for 
several x values and compared with the standard. We can see that the difference 
between the standard and a; = 0.1 is of order 10~* or less, but it becomes even 
more important while x increases. 

In Figure [531 we plot the time and temperature evolution of the mass fraction 
of several light elements plus metals (by metals we mean the sum of elements 
heavier than ^He, that is ^Li, Li and so on) in the ordinary case and in 
both sectors in the mirror model; we used for these calculations the parameters 
X = 0.7 and (3 = 5. We can see that BBN in the mirror sector is much more 
different from the standard than the ordinary sector one. This is a consequence 
of the higher number of DOFs. 

In Figure ESI we plot the final abundances of '^He, D, ^He, '^Li and metals 
in the ordinary sector as a function of x. We can easily infer that for x < 0.3 
the light element abundances do not change more than a few percent. 
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Figure 5.6: Light elements produced in the ordinary sector during 
BBN. On every plot we report the elements {^He, D, ^He, Li and 
metals) for several x values, compared with the standard and the 
"order 0" quadratic approximation of AA^,^ in equation (IS.lSp . 
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5.4 The light elements production in the mir- 
ror sector 

Even mirror baryons undergo nucleosynthesis via the same physical processe 
than the ordinary ones. But the impact of the ordinary world on the mirror 
BBN is dramatic since the ordinary contribution on the number of total mirror 
DOFs scales as ~ - see equation (I4.23p . The knowledge of the element 
abundances in the mirror world can be important for the study of its evoultion 
in later epochs as the galaxy formation, and in particular, for the formation and 
evolution of the mirror stars [82]. 

We can use the Kawano code also for the mirror nucleosynthesis; once again 
there are two different ways to input the neutrino number N,y: the order zero 
AA*",^ = const oc and the numerical values obtained solving the equations 
with Mathematica. The physical parameters t], t and h have the values in 
equation fl5.18p . 

An other parameter which affects mirror BBN is the mirror baryon density 
{fl'^ ~ 1 ^ 5Qb), which raises the baryon to photon ratio rj' - see §2.5j in our 
program Q'^ is introduced in terms of the ratio P = Q'^/Qb - see equation (I4.33p . 

In Table 1^751 the program output, that is the light elements produced in the 
mirror sector at the end of BBN process (at T ~ 8 ■ 10~^ MeV), is reported for 
several x and (3 values and compared with the standard. 

In Figure [5?7I we plot the final abundances of '^He, D, ^He, '^Li and metals 
in the mirror sector as a function of x. Since T' ~ xT, at higher x the mirror 
sector becomes hotter; hence, we expect to have mirror abundances closer to 
the ordinary ones at higher x. Moreover, we note that in general lower j3 implies 
final mass fractions closer to the standard. 

A very rough approximation for the mirror '^He production can be achieved 
from the following argument: for any given temperature T', using equations 
f l4.23p and 04.241) and assuming a ~ 1 we have 




i/(r)^ 5.5^1 + -J— (5,19) 

for the Hubble expansion rate. Therefore, comparing H{T') with the reaction 
rate r(T') oc T'^ - see equation (12.8 ip - we find a freeze-out temperature 

Tly = {l + x-y/'Tw (5.20) 

which is larger than Tw, whereas the time scales as 

t'w = tw/il + x-^f^ <tw (5.21) 
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Table 5.5: Light elements produced in the mirror sector. 



elements 


X = 0.1 {13 = 5) 


X = 0.1 (/3 = 1) 


X = 0.3 (/? = 5) 


n/H 


5.762 -10-25 


8.888 -10-1^ 


2.590 -10-22 


P 


0.1735 


0.1772 


0.3646 


D/H 


1.003 -10-12 


1.331 -10-6 


4.838 -10-9 


T/H 


9.679 -10-21 


3.068 -10-'^ 


1.238 -10-13 


^He/H 


3.282 -10-6 


5.228 -10-6 


3.740 -10-6 


^He 


0.8051 


0.8226 


0.6351 


^Li/H 


7.478 -10-21 


8.638 -10-15 


1.309 -10-1^ 


^Li/H 


1.996 -10-^ 


5.712 -10-s 


3.720 -10-* 


^Be/H 


1.996 -10-^ 


5.711 -10-*^ 


3.720 -10-* 




^.oo^ -iu 


z.uou -iu 


c, nop. . 1 n^ii 
o.yzu -iu 


elements 


X = 0.3 {(3 = 1) 


X = 0.7 {/3 = 5) 


X = 0.7 {(3 = 1) 


n/H 


1.915 -10-16 


1.726 -10-19 


2.076 -10-16 


P 


0.3675 


0.5924 


0.6017 


D/H 


7.094 -10-6 


3.279 -10-' 


2.235 -10-5 


T/H 


2.190 -10-^ 


3.722 -10-1° 


7.328 -10-* 


^He/H 


6.880 -10-6 


4.691 -10-6 


9.719 -10-6 


^He 


0.6326 


0.4077 


0.3984 


^Li/H 


1.660 -10-1^ 


3.361 -10-16 


1.951 -10-1^ 


'^Li/H 


8.930 -10-9 


7.962 -10-9 


1.120 -10-9 


^Be/H 


8.878 -10-9 


7.962 -10-9 


1.064 -10-9 


^Li + /H 


2.514 -10-12 


3.949 -10-13 


1.814 -10-1^ 


elements 


x = 0.5(/3 = 5) 


X = 0.5 {p = 1) 




n/H 


1.840 -10-20 


2.058 -10-16 




P 


0.4966 


0.5028 




D/H 


6.587 -10-s 


1.352 -10-5 




T/H 


2.108 -10-11 


4.358 -10-* 




^He/H 


4.172 -10-6 


8.232 -10-6 




^He 


0.5035 


0.4974 




^Li/H 


1.016 -10-16 


1.790 -10-1^ 




'Li/H 


1.535 -10-« 


2.948 -10-9 




^Be/H 


1.535 -10-s 


2.891 -10-9 




^Li + /H 


3.827 -10-12 


1.657 -10-13 
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Figure 5.7: Light elements produced in the mirror sector during 
BBN. On every plot we report the elements {'^He, D, ^He, Li and 
metals) for several x values, compared with the ordinary sector stan- 
dard and the "order 0" quadratic approximation of in equation 
( IS.lSp . For every x we used two different values for the barion ratio 
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(obtained using the relation t oc H^^). In addition, rj' is different from 77 ~ 
5 X 10^^°. However, since T/v depends on baryon density only logarithmically 
(see eq. (12.751) ). the temperature remains essentially the same as T^r, while 
the time t'^ scales as 

4 = W(l + x-^)'/' (5.22) 

Thus, for the mirror ^He mass fraction we obtain - see equation (12.871) : 

2 exp~*^/^(i+''''')''" 
- 2X(4) = X ^ Q^-pQ/T^{i+x~^Y/^ (^-23) 

We see that is an increasing function of In particular, for x ^ one has 
1^4 1. Hence, in this case is always bigger than Kj. In other words, if dark 
matter of the Universe is represented by the baryons of the mirror sector, it 
should contain considerably bigger fraction of primordial ^He than the ordinary 
world. In particular, the helium fraction of mirror matter is comprised between 
20% and 80%, depending on the values of x and i]' . This is a very interesting 
feature, because it means that mirror sector can be a helium dominated world, 
with important consequences on star formation and evolution, and other related 
astrophysical aspects. 



Chapter 6 

BBN when the mirror sector is 
millicharged 



If the mirror particles are milli- charged with respect to the ordinary photon, 
see chapter |H there can be processes, hke the pair annihilation and produc- 
tion, which lead to entropy and energy exchanges between ordinary and mirror 
sectors. 

A consequence of this is that the entropy densities ratio x we defined in §4.31 
is not a constant of time anymore, but changes during the universe evolution. 

In this chapter we analyze these processes and work out rough bounds on 
the photon mixing parameter e; finally we will critically discuss them in order 
to propose further developments of the work began in this thesis. 

6.1 Physical processes involving mirror and or- 
dinary matter 

The processes of lower order involving mirror and ordinary matter are propor- 
tional to e^. There are three: the pair annihilation-production e~^e~ ^ e'^e'~ , 
the scattering ee' ^ ee' and the plasmon decay 7 e'^e'~ . It can be shown that 
the plasma process is negligible: the plasma frequency in a relativistic plasma 
indeed is [94]: 



Since we are interested in particles having the same or higher mass than the 
electron, the plasma process is below its threshold (2me) and can be ignored. 
The other two processes are analyzed below. 




(6.1) 



87 



§6. J. Jl Pair annihilation and production 



88 



6.1.1 Pair annihilation and production 

The cross section for the process e^e^ ^ e'^e'^ is worked out in detail in 
appendix IB.2I The total cross section for the special case rrie = rrie' = m is: 

a = -Tia^ — (s + 2m^) 6.2 
3 s'^ ^ ^ 

where the calculation has been performed in the center of mass frame of reference 
and s stands for the total four-impulse squared: s = \pi - see IA.5I 

Once the cross section is known, F can be worked out; at low temperature 
we can follow Gondolo and Gelmini [90], who demonstrated that the thermal 
average can be achieved solving a single integral: 



2 ^m*TK^[ 



X 



where a; = ^, i^'j are the modified Bessel functions of order i and the factor 
I arises because the initial particles are not identical. This formula is valid 
for particles with Maxwell-Boltzmann statistics or having any other statistics 
provided that T < 3m. The corresponding quantity in the center of mass (CM) 
frame is: 



{(Tv) 



CM 



21 NT-l 



{(yvM<f,i) (6.4) 



Kl{x) 

Recalling n from equation (12.331) . F is finally obtained as: 

F = {avf^'n (6.5) 

At high temperature we can use an ultra-relativistic approximation, which is 
described below. 

Ultrarelativistic approximation 

For ultrarelativistic (UR) fermions - see equation flA.22p - we have: 

4 2 rr.2 

(^UR = — s ocT 

3 s 



n{T)uR = ^^J' v = c = l (6.6) 

4 TT^ 



This way we get the ultrarelativistic interaction rate: 



= n{ca) = ^e'a'T ^ 0.2e'a'T (6.7) 
III 
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Figure 6.1: F fore = 10^^ calculated from Gondolo and Gelmini for- 
mula Ii6. 5\) (red) and in the ultrarelativistic approximation in ( [g. ?] ) 
(blue). F is in units log(F/MeV) 



In Figure 16.11 F is plotted for e = 10~^ in comparison with the ultrarelativistic 
approximation. We can see that the relative difference between the two formulae 
is of order 10% or less in the range 1 < T < 5MeV, but it becomes much higher 
for any other temperatures. 

Comparing the interaction rate (16. 7p with the Hubble parameter H - see 
equations (12.511) and (12.461) - we can work out a rough bound on the millicharge 
e. Indeed, to preserve mirror matter out of equilibrium when the neutron-proton 
reactions (I2.8UI) freeze-out, that is at T = ~ 0.8MeV^ - see equation (12.811) 
- we must have 

— = [ea] —— < 1 6.8 

H 27rl.667^^ ' Tw 

from which we obtain the bound on the photon kinetic mixing parameter 



e< 571 J— -4.7-10^9 (6.9) 

- V Mp ^ ' 
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6.1.2 Scattering 

The calculations we made to calculate the cross section for the process ee' ee' 
are reported in appendix IB. 31 The differential cross section is: 



da 

dn 



2E^ +p^[l + cos^e)] + AE^p^ cos' (^^^ - 2m\E^ - p^cosO) + 2m^ 

e2«2 1 ^ ^ 

(6.10) 



32sin4(f) E^p^ 



Clearly, the differential cross section and its integral diverge at small angles; 
anyway this process can be ignored at first order in x, as we can see from the 
following argument. 

Let us consider for simplicity relativistic particles, so that n oc T^, T oc T 
and (E) oc T; to compare the scattering process with the pair annihilation, we 
can use the order- zero approximation for the energy loss - see equation fl6.16l) : 

~ -rne±(^e±) ocT^ (6.11) 

dt 

For the process ee' ^ ee' instead we have s oc TT' from which follows F oc 
rigis oc YT and therefore 

^ ^ -Vn,,,{E,-,) oc ^r'T' ^ '^x' (6.12) 

BBN bounds impose x < 0.64 and therefore we can ignore at first order 
since 

dpscat ^ ^ dp 

ann -, o\ 

6.2 Energy transfer between the two sectors 

In this section we will work out the energy exchange between the mirror and 
the ordinary sectors using various approximations. 

For any processes the number of interacting particles N changes in agreement 
with 

f = ("^' 

where F = {o-v)ntarget as usual and N = nV, oc a^. 
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6.2.1 Neglecting the universe expansion 

The roughest approximation we can make is to neglect the universe expansion, 
that is to assume a = 0; in this case we will have: 
dn ^ dp dn 

If ordinary particles, which are hotter that the corresponding mirror ones, an- 
nihilate by this process creating mirror matter and viceversa, the change in the 
ordinary world energy density will roughly be: 

dpo 
dt 

= -T{T)no{T){Eo{T))+T{T')nUr){Em{r)) 



dno_^{T) 
dt 



{Eo{T)) 



drim^oiT') 
dt 



{Em{T')) 



- - r(T)p„(T) + -r(T')p„(T') 



TT 

30 
30 



g{T)T{T) 



g'{T')V{T') 



g{T)T{T)-g'{T')T{T') 



(6.16) 



where the index o stands for "ordinary" and the index m stands for "mirror". 
In the previous calculation the thermodynamic functions - p, (E), g - are not 
of the whole universe but just of the involved particle species. 

We stress that the term describing the energy transfer from the mirror to the 
ordinary sector is of order ~ < 0.15; thus, considering the present accuracy 
level, we can neglect it . In the following we will consider the only energy transfer 
from the ordinary to the mirror sector. 



Order of magnitude of the energy loss 



Using the ultrarelativistic interaction rate for the pair annihilation in equation 
( 16.71) in the order-zero approximation for the energy density loss (I6.16p we get 



dp 
'dt 



X 



-3 ■ l{}~^e^g{T)T^ 



\ — X 



9{T) 
5 9'iT') 
9iT) 



- -3 ■ l{]-^e^g{T)T^ (6.17) 

The substitutions IMeV ~ 1.6 ■ lO^^s"^ and g{T) ~ 4 will finally lead us to 

dp 



dt 



2 ■ lO^^e^T^ {MeV s)'^ 



(6.18) 
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which is the same energy loss used by Carlson and Glashow [55]. We can use 
this formula to work out a bound on e more accurated than the one in equation 
(16 ■91) . Indeed, the energy transferred to the mirror sector is about: 

Ap'(T) = / ^dt' = 0.602 g,± ■ lO'^e^^-^/^ ^^r^s ^q ^q^ 

J (JjIj 

where we used the substitution in equation (12.471) and ignored the g dependence 
on T, that is 

t ~ 0.301 g-'^^ ^ dt ^ -0.602 g-'^' ^ dT (6.20) 

We recall also that Ag = 2 ■ ^ ■ AN^ = 1.75 ■ AA^^ and thus 

f^<H5^^0.163A7V„ (6.21) 
p 10.75 

But, since p' > Ap' , we have 

^ = ^ ^ 0.6 ■ 10-^^ < ^ < 0.163 AK (6.22) 
p gy/g T p 

In conclusion, e is bounded to 



e< jl^^^lO^^-V^ (6-23) 
- V 6 30 ge± 0.602 Mp ^ ^ ^ 

where as usual Tw ^ 0.8 MeV, Mp ~ 1.2 ■ 10^^ MeV; at Tw we have - see gSXl 

g ~ 10.45 10.45 <g< 13.0 ge± ~ 3.39 (6.24) 

Using these values and assuming the conservative bound AN,j = 0.5 we finally 
get the bound on the kinetic mixing parameter: 



e < 0.246 ■ 10"^ ^y^AiV^ < 0.80 ■ lO^^AiV^ < 5.7 • 10"^ (6.25) 

6.2.2 The number density in the expanding universe 

Let us now take in account of the well-known fact that the universe expands 
and consider the ordinary particles annihilation followed by the mirror particles 
creation, that is the process e~^e^ ^ e'~^ e'~; we inferred in §6.11 that this is the 
most important process involving particles from both sectors. 
We can use equation (I6.14p . which still holds and implies: 

^ = -TN, and ^ = TiV, (6.26) 
dt dt 
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Making the substitutions 

Ne' oc rie'a^ Ne oc n^a^ H =- (6.27) 

a 

we get the Boltzmann equation 

he' + Sue'H = Trie (6.28) 
Let us now rewrite equations (12.461) and (I2.47P as: 

t - a—^ H ~ b (6.29) 

r2 Mp ^ ^ 

where a and b are constants and we neglet the g dependence on T. Equation 
(I6.28P becomes in terms of the temperature: 

dup' Up' 

^-3^ = /(r) (6.30) 

where 

, -2aT(T)nJT) Mp 

f{T) = ^ !^;^ ' (6.31) 

The first step to resolve the differential equation (I6.30p is to find a solution to 
the homogeneous associated equation, which gives: 

n^, = cT^ (6.32) 

where h stands for homogeneous. Then we make the constant c vary in order 
to find the general solution. Assuming c = c(T) in (16.320 and substituting it in 
(I6.30p we finally get 



c(T) =c,+ r ^dy (6.33) 



where ci is a constant. Therefore the general solution of (I6.30p is 



ne'{T) 



ci+ / -^dy 



(6.34) 



Initial condition: empty mirror sector At high temperatures - such that 
T 3> me - electrons are ultrarelativistic and we have rie oc T^, F oc T - see 
§6.1. It it follows that /(T) = C2T where 

3 ■ 0.301 (^1 (yfei 2 2 r r i n-7 9e^ , _^ 2 /ro^\ 

C2 = 5 — Mpae~ 5.6 -10 ■—=Mpe (6.35) 
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Solving the integral in (16.341) we finally get 

ne'{T) = ciT^ + C2T'^ (6.36) 

Let us impose now the initial condition requiring that the number of mirror 
partilces can be neglected respect to the number of ordinary particles, that is, 
the mirror sector is empty: 

hm — — oc hm -—- = 6.37 

T^oo ne{l ) T-^co i i 

This implies ci = and therefore the general solution of equation (I6.30p becomes 

nAT) = T^ r^dy (6.38) 
which becomes, in the ultrarelativistic limit 

He' = CsT^ (6.39) 

Initial condition: x dependence An other possible initial condition is that 
at high temperatures the mirror fermions / are at equilibrium with the other 
mirror species, that is: 

lim rie' = - -44T'3 that is lim — = — (6.40) 

T^oo 4 TT-^ T^oo Tie T'^ 



This implies that 



4 vr^ 4 TT^ 



from which follows 

c. = \%^x' (6.42) 



6.2.3 The energy density 

The energy density of mirror particles can be written approximately as 

Pe'^ne'(^e') (6.43) 

while the energy transfer from the ordinary to the mirror sector is 

Ape' = Ar2e'(^e') (6.44) 

We worked out ng' as a function of the ordinary sector temperature T in §6.2.21 

In the ultrarelativistic limit we have Ang' = C2T^, (-Eg') — 3.15T and we can 
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use these formulae to work out the energy transfer from the ordinary to the 
mirror sector: 

Ape'(T) ^ 3.15 caT^ (6.45) 
From which follows, applying the same procedure of §6.2.11 
An' 1.8 ■ 10-6 ■ ^ Mp e^T^ n' i 

^ ~ , ^ < ^ < -AN, (6.46) 

Using again fl6.24p . we get the bound on the kinetic mixing parameter: 



'1 T^'^gs/g 1 ,^.T^ 



e<J- —IXJL ^ 106^ . Jan, - 3.3 ■ 10-9 (6.47) 
- V 6 30 ge± 1.8 Mp ^ ^ ' 

6.3 Prospects for future researches 

In the previous chapters of this thesis we analyzed the features of the mirror 
model and worked out in detail its primordial nucleosynthesis if millicharged 
interactions are not present. 

Future researches should focus on BBN in presence of these interactions to 
work out more precise bounds on e and x. In particular, the processes we an- 
alyzed in §6.11 lead to energy and entropy exchanges between the sectors and 
therefore to a not trivial evolution of the entropy ratio s' / s with time. In gen- 
eral we expect millicharged interactions to raise the mirror sector temperature, 
energy and entropy and lower the corresponding quantities in the ordinary sec- 
tor, making s'/ s grow with time till the temperature becomes smaller than the 
threshold of the involved processes, that is T ~ me/2. 

It is also important to stress that the presence of the mirror universe and 
its being a candidate for dark matter does not affect only BBN, but also other 
cosmological and astrophysical processes. The parameters x and (3, which are 
enough to unambigously fix the symmetric model, must be therefore chosen such 
that all these processes take place and are compatible with the experimental 
observations. 

For instance, the structure formation in presence of mirror dark matter has 
been analyzed by Berezhiani et al. and the previsions of this model on the 
cosmic microwave background (CMB) and the large scale structure (LSS) have 
been compared with those of the cold dark matter (CDM) model [81,83]. 

What came out from these comparisons is that, if x < 0.2, CMB and LSS 
power spectra are equivalent for the CMB and CDM models and it is possible 
that all DM is made of mirror baryons, while for higher x these spectra strongly 
depend on the amount of mirror baryons, which must, to save the compatibility 
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with data, constitute not more than 20% of DM. This happen because decou- 
phng of matter and radiation must take place before in the mirror sector and 
be complete before the ordinary one begins. 

Let us now recall the three rough bounds on e we worked out in equations 
( 16.91) . (16.251) and (16.471) . All them have the same order of magnitude, that is 
e < (3.3 -r- 5.7) X 10^^ and are proportional to y/AN^, §. We recall also from 
equation (I4.27P that AN^ oc and hence 



Therefore we can see that imposing x < 0.2 implies AA^^, < 10 ^ and hence the 
bounds on e must be scaled by a factor 



and hence they become about 7 times smaller, so that e < (4.6 ^ 8.1) ■ 10^^'^ 
These new bounds are not compatible with e ~ 5 ■ 10~^ required to explain the 
DAMA signal in terms of mirror baryons scattering - see §4.51 

Nevertheless, there are at least two ways out: the first is found assuming 
that DM is not only done by mirror baryons but also by some other kind of cold 
matter; this way x is only bounded to be less than 0.64 by BBN and we have 
the bounds in equations (16.91) . (I6.25P and (16.471) . 

The second is based on the asymmetric mirror model introduced in §4.4[ 
where m'^ ^ 5mB, rrie' ~ 100 mg. In this scenario mirror matter and radiation 
decouple before, thus the only bound on x arises from BBN, that is x < 0.64. 
Moreover, since the threshold for the process e~^e^ ^ e'+e'^ is at Tfhr ~ ^ ~ 
25 MeV, the energy and entropy exchanges between the sectors lose their efficacy 
before. 

In particular, at T < Tt^r exchanges between the sectors are not efficient 
anymore and hence the universe evolution goes on just as in chapter [5l the ratio 
s'/s is a constant, both x and AN^, are approximately constant (if we neglet 
the small variation due to the a factor in equation (I4.25P ) and so on. Therefore 
we can replace Tw by Tthr in the three bounds on e we worked out, that is 
equations (16. 9p , (I6.25P and (16.471) , and this way they become higher by a factor 




(6.48) 



V2 ■ 10-2 ~ 0.14 ~ - 



(6.49) 




(6.50) 



that is, we have 



e < 2^3.5 



10 



,-8 



(6.51) 



Cosmological limit on e will be discussed in more details elsewhere [95] 
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In such an asymmetric mirror scenario we may have enough baryons to consti- 
tute all the dark matter required by experimental observations of the universe. 
At the same time, this scenario is compatible with CMB and LSS and its BBN 
impose bounds on the kinetic mixing parameter e of order a few unity 10~®. 
This asymmetric model may also explain the DAMA/Nal annual modulation in 
terms of scattering of millicharged mirror baryons on ordinary matter, since e is 
high enough. But all these intuitive hints should be developed in a more rigor- 
ous way to work out the parameter evolution with time and the e value required 
to explain the DAMA/Nal modulation and finally verify the compatibility of 
the model with all the experiments. 



Conclusion 
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It is nowadays well known that particle physics and astrophysics are tighly 
bound; indeed, a good knowledge of the interactions between particles can help 
in understanding or building models for cosmology and astrophysics, as well 
as the present and primordial universe provide high energy environments which 
are unreachable in laboratory experiments but are crucial for studies concerning 
elementary particle interactions. This is why in the recent years a new branch 
- called astroparticle physics - has became even more popular. 

This thesis is devoted to some astroparticle implications of a mirror sector, 
which is a parallel hidden sector of particles and interactions, completely iden- 
tical to the ordinary (observable) particle sector, i.e. standard model and its 
possible extensions, and coupled to the latter via gravity. In other words, the 
mirror sector has identical particle physics to that of the ordinary one. The 
mirror sector presence restores the apparently unexplicable lack of the left-right 
symmetry in nature; moreover, mirror baryons are a viable dark matter can- 
didate, since they couple to the ordinary matter mainly via gravity and their 
density can be as high as the dark matter one should be (about 5 times bigger 
then the ordinary one). Therefore experiments searching for dark matter (such 
as DAMA) are natural candidates for detecting mirror matter. 

The presence of the mirror sectore does not introduce any new parameters 
in particle physics (as soon as mirror sector particles and interactions are iden- 
tical to the ordinary ones). However in cosmology they should have different 
realizations; in particular the temperature of the mirror sector T' should be less 
than the ordinary one T, otherwise it would spoil standard BBN predictions, 
which are one of the key stones of the standard cosmology. Therefore the value 
X ~ T'/T is a free parameter which value is restricted by BBN conditions to 
X < 0.64. 

An interesting possibility is that any ordinary particles may be mixed with 
its mirror twin. In particular ordinary photons can have kinetic mixings with 
mirror photons, as well as ordinary neutrinos, ir^, neutron and so on can be 
kinetically mixed with their mirror partners. 

This thesis is devoted mainly to the kinetic mixing of ordinary and mirror 
photons and its cosmological consequences regarding BBN epoch of the primor- 
dial production of light elements and is constructed in the following way. 

The first two sections are devoted to basic concepts and models of particle 
physics and cosmology. We start with a brief review of the standard model 
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of elementary particles and their interactions and the standard cosmological 
paradigm, followed by a discussion about the origin of the electric charge quan- 
tization in particle physics (third chapter), which is induced in the quantum 
standard model by the anomaly cancellation conditions and naturally arises in 
grand unified theories or in presence of magnetic monopoles. 

In particular, it is a challenging question whether " milli" charged particles, 
having very tiny electric charges could exist in the universe. This possibility is 
present in models where two photons kinetically mixed are present, what can 
be naturally achieved in models containing a mirror sector, in which case mirror 
electron, positron, proton and so on would emerge as millicharged with respect 
to the ordinary photon. Of course, if mirror particles are millicharged, this 
would open new possibilities for their detection. 

In chapter [5] we worked out a detailed analysis of BBN in the mirror model 
as well as in the mirror sector itself. To perform these calculations we first found 
a system of equations which numerical solution gives the mirror temperature T' 
corresponding to any ordinary temperatures T; once these temperatures were 
known, we were able to work out the number of degrees of freedom and the num- 
ber of extra-neutrinos in both sectors. A special feature of models containing 
the mirror sector is that the number of equivalent neutrinos (that is the total 
number of degrees of freedom expressed in terms of neutrinos) changes during 
the universe evolution. Finally we modified the Kawano code for BBN to let it 
read the number of neutrinos as an external input changing with time. These 
calculations had never been done before with the same accuracy. 

In chapter [H] we insterted the kinetic mixing of ordinary and mirror photons 
and analyzed the main processes leading to energy and entropy exchanges be- 
tween ordinary and mirror sectors. We made some extimates of the maximum 
value for the photon mixing parameter e, which came out to be a few units 

An interesting possibility is that mirror baryons scattering on ordinary baryons 
may explain the DAMA/Nal annual modulation if e ~ 5 ■ 10~^, that is of order 
of our upper bound. 

We would like now to stress some final remarks which constitute a starting 
point for further future researches. In our extimates the maximum value of e 
(let us call it emax) is proportional to x^; the structure formation bounds x to 
be smaller than 0.2 and this leads emax to be about one order of magnitude 
stronger, that is, not compatible with the DAMA/Nal signal. 

Anyway, the bound x < 0.2 applies only if all dark matter is made of mirror 
baryons; hence we can assume that at least 80% of it is made of some other 
kind of cold matter and save compatibility with DAMA. 
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Otherwise we may use an asymmetric mirror model, where m'^ ~ Smg, 
m'g ~ 100 mg. In this scenario the condition Q!^ ~ SQ^ is naturally achieved if 
mirror and ordinary baryons have the same density. Moreover, having heavier 
mirror electrons lead to an earlier freezing of the pair annihilation-production 
processes and this allows Cmax about one order of magnitude bigger (~ 3- 10~^), 
that is, compatible with the DAMA interpretation. 

In conclusion, the mirror model is a challenging theory, supported by theo- 
retical symmetry reasons as well as by its providing a natural candidate for dark 
matter. The mirror sector interacts with the ordinary one via gravity, but it 
may also have small interactions via photon mixing (the mirror model also pro- 
vides a natural explication for this smallness). Observable effects of the mirror 
matter may be found in experiments searching for dark matter as well as in the 
primordial universe. In particular, calculations of the light elements produced 
during BBN should be performed with higher accuracy for both symmetric and 
asymmetric mirror models and the compatibility of the parameters with the 
DAMA/Nal annual modulation should be verified. 
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Appendix A 

Useful formulae for cosmology 
and thermodynamics 



A.l General Relativity 

The essence of Einstein's theory is to transform gravitation from being a force 
to being a property of space-time, which may be curved. The interval between 
two events 

ds"^ = g^ij{x) dx'^dx'^ (A.l) 



is fixed by the metric tensor g^p which describes the space-time geometrjt 
{g^^Qf^jS = 5^). For the Riemannian spaces, the tensor of curvature is 

where the F's are ChristofFel symbols 

1 



'^P i^a I p/i per pM po- (' A O^l 



J- a/3 — 2 ^ 



dx^ dx"" dx'^ 



(A.3) 



The equation of motion of a free particle is determined by the space-time metric 

d-x- d^<W^^ (A.4) 



ds ds 

SO that free particle moves on a geodesic. 

On the other hand, the metric gajs is itself determined by the distribution 
of matter, described by the energy- momentum tensor T^p, according to the 
Einstein equations 

Gap = Rai3 - ^Rgap = SnGT^fB + Ag^p (A. 5) 



"'^ Repeated indices imply summation and a, j3 run from to 3; x'^ = t is the time coordinate 
and {i = 1, 2, 3) are space coordinates. 
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gAj Flat models 
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where 



Ra 



:/3 



:/3 



(A.6) 



are respectively the Ricci tensor and the Ricci scalar, and A is the cosmological 
constant. For the FRW metric, the non zero components of the Ricci tensor 
and the value of the Ricci scalar are 

a2 



Rqo 
R = 



-3^ 


Rij 




a 


6 


a a? 


k' 


- + - 


+ - 




a 





+ 2— + 



2k 



9ij 



For a perfect fluid the energy-momentum tensor has the form 

Tal3 = {P + p)UaUfS - PQafi 



(A.7) 



(A.8) 



where p and p are respectively the energy density and pressure of the fluid, 
Ua = QajS dx^ /ds IS the fluid four- velocity. Considering the symmetries of the 
FRW metric (uniformity and isotropy), which demand that = 1 and = 
in the comoving coordinate system, we obtain Tap = diag{p, —p, —p, —p). This 
is valid also in presence of the cosmological constant, if we substitute p and p 
as indicated in the following 



p 



P-PA , 



p ^ p + pA] 



A 



PA 



-PA 



8nG 



(A.9) 



Therefore, the Einstein equations for a Universe described by the FRW metric 
are reduced to equations fl2.10p and (12.121) . where the relations between the 
energy and pressure densities are in general related as p = wp. In particular, 
for the dominance of relativistic and non relativistic matter we have respectively 
w = 1/3 and w = 0, while for the vacuum energy dominance one has w = —1. 

A. 2 Flat models 



Here we report two special cases of the relationships (12. 221) - (12. 271) : dust or 
matter dominated Universe [w = 0) 



a{t) 

t 
H 



|-i/o(l + .) 

1 

2 



3/2 



(A.IO) 

(A.ll) 
(A.12) 

(A.13) 



§A.3I Thermodynamic functions in special cases 
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tn 



5Ho 
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6nGf^ 

and radiation dominated Universe {w = 1/3) 



a{t) 

t 
H 
% 
to 

Pr 



ao 
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to(l + ^) 

l = Hoil + z) 



2Ho 
3 



327rGt2 
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(A. 14) 
(A.15) 

(A.16) 

(A.17) 
(A.18) 
(A.19) 
(A.20) 

(A.21) 



A. 3 Thermodynamic functions in special cases 

In this section the number density n, the energy density p and the pressure p 
functions in some special simpler limits are reported. 



Ultrarelativistic limit: T ^ m , T ^ fi 



p = 3p 



8 - ^ 



■/30 QspirijiT'^ 
■/30 QspiriA^ 



(BOSE) 
(FERMI) 



n 



C(3)/vr2^.p.„^T3 



(BOSE) 
(FERMI) 



(A.22) 



Non-relativistic limit: T <^ m, 



n 

P 

P 



g spin A I ^ ) ^ 

-(m-M) 

mn oc exp t 



ni = <^ p 

m 



(A.23) 



A. 4 Entropy: definition and conservation 

The second law of thermodynamics states that 

TdS = d{pV) +pdV -p d{nV) = d[{p + p)V] - Vdp - p d{nV) (A.24) 
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where p and p are the equihbrium energy density and pressure; differentiating 
equation fl2.35p with respect to T we get 

dp g 1 



dT Gn^Jo E (expVil^^ 



exp T 



E d / n 
T^^dTVT 



dk 



where we used 



9 



Gvr^ Jo dk 



'E d /fi 



dk 



f 



exp T ±1 



(A.25) 



dl 
dk 



ET 



2 _ E- p 

J rj-t 



E"^ = k'^ + 



(A.26) 
(A.27) 



dE _k 
~dk ~E 

Integrating equation (1A.25P by parts and using equations (I2.34p and (12.351) we 
get: 

^ k=oo 



dp 
dT 



9 



67r2 



k'T-f- — + 



E d / n 



T2 dT \T 



k=0 



+ 



_9_ 

Gtt^ Jo 



P+P , d f P 



Using this equation in flA.241) we get: 

P + P 



(A.28) 



dS 



^d[{p + p)V]-^V 



— + nT.^(^ 
T dT \T 



dT - 



jidijiV) 
T 



d 



V ■ 



p + p — np 



T 



(A.29) 

that is equation (12.361) if we assume as usual V (x B? . In the following we will 
always assume p <^ T. To demonstrate that S = constant we recall the first 
law of thermodynamics (12.171) : 

d{pV) = -pdV - pdinV) (A.30) 

Sobstituting this identity in (1A.29|) and using again (1A.28P we finally get: 

dS = Q (A.31) 
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A. 5 The four-impulse average for massless par- 
ticles 

Let us consider a process where two particle 1 and 2 with four impulses Pi = 
{Ei, Pi), i = 1,2 interact. When they are ultrarelativistic the interaction rate 
for processes having a oc s - see §B.1I - is 

r = n{va) ~ n{a) oc n{s) (A. 32) 

It is therefore useful to work out the four-impulse squared average; the total 
four- impulse squared is defined as s = |pi +p2p and its average on the impulses 
is 

/ V ^ / |Pl +P2|Vl(pi,t)/2(p2,t) d^Pl d^P2 . . . 

Jfl{pi,t)Up2,t)d^Pld^P2 ^ ■ ' 

where fi{pi,t) d^pi has been defined in equation fl2.3ip . For massless particles 
we have: 

Pi = (|Pi|,Pi) ^ bi+P2r = 2|pi||p2|(l-cos^) (A.34) 

/, oc ^ (A.35) 

exp T ±1 

j d^pi d^P2 = 4:71 J |pi|2rf|pi| ■ 27r|p2|^y" rf|p2| j^^ dcosO (A.36) 
The denominator of flA.33P is: 



/ /i(Pi,t)/2(P2,t) rf^Pi d'p2 = {Anf f / (A.37) 

exp T ±1 exp T ±1 

and can be easily worked out recalling that 

f Pld\pi\ 

/ |p-| '^i ultrarelativistic (^/i.Ooj 

exp^ ±1 



where ultrarelativistic IS givcu in (]A.22p . Therefore equation (1A.37P reduces to 



^4 2^2^6 y / 1 (for any present BOSONS) , . . 

D47rc3i (for any present FERMIONS) ^ 

Similarly the numerator of flA.331) turns out to be proportional to the energy 
density p for ultrarelativistic particles. Using flA.34p we have: 

167r2 r\l - cos9)dcos9 [ 1| d|pi| / 1| ci|p2| (A.40) 

•^"^ exp T ±1 exp T ±2 
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We can observe that 



Pi 



'^|Pi| Pi ultrarelativistic 

(A.41) 

exp T ±1 

where pi uitrareiativistic IS given in (]A.22p . Finally equation (1A.40P becomes 



32 



225 

Therefore (s) finally is: 



450CI 



X 



Hence 



14.6 



(for any present BOSON) 
(for any present FERMION) 



(for any present BOSON) 
(for any present FERMION) 



(boson — boson scattering) 



(A.42) 



(A.43) 



17.1 



(boson — fermion scattering) 



= 19.9 



(fermion — fermion scattering) (A. 44) 



Appendix B 
Cross sections 



In many books, such as [31,32], we can find expressions for the differential cross 
section Let us consider a scattering process in which two particles - leptons 
or photons - with four momenta pi = {Ei, pi), i = 1,2 collide and produce 2 
final particles with moments p'j = {Ej,Pf), f = 1,2. In the center of mass 
frame (CM) we have = —p'2 and the differential cross section is 




nK2mO Ip^i 

647r2(Ei + E2)2|pi| 



X (B.i; 



where / runs over all the present lepton species, having masses mi and X is the 
spin average: 

X = lj:\M\' (B.2) 

^ r,s=l 

Below we will work out ^ for some processes of interest for this thesis. 

B.I Neutrino scattering 

In this section we will assume fc^ -C Tn"^, k"^ <^ m| where k is the intermediate 
boson momentum, so that the weak boson propagators reduce to the Fermi 
constant. 

Let us consider the neutrino scattering z/^a z/^a with a = e, fi,T mediated 
by the charged weak bosons as depicted Figure IB. II A. We will indicate 
below by the apex CC when referring to charged currents. The cross sections 
we will get are the same we would have for the scattering Vafl ^ ^a(i- The 
Feynman amplitude for the processes we are considering is: 

\Mfc ^ ^ [u{k')^^{l - 75)«(p)] W)i,{l - i,)u{k)] (B.3) 
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from which follows 

X^^ = ^tr [7^1 -75)(j^ + m)7^(l- 75)^'] tr [7,(1- 75)^7.(1 -75^'] = 

= ^256 {p,kn{plk'n (B.4) 

where Gp — 1-1664 ■ lO^^GeV^^ is the Fermi constant. The differential cross 
section for this process is therefore 



dVtj ~ Un^s \p\ 2 ^""^ 



^ 128 ip'X) (B.5) 



There are two extreme cases of interest in physics: the first is an interaction 
of a neutrino with an ultrarelativistic target. In this limit Ei, ^ Ma] therefore 
s = 2pf^k^ = 2p'^k'^^, p' = p and 

/ do'\'^^' G% nn da G"^ 



The second is an interaction with an ultramassive target, when Ey <^ M^, p ^ 0; 
therefore 

s ~ Ml p^k^ = p'^k'^ = MaE, (B.7) 

Let us now turn to the neutrino scattering ^ with i,a = e, fi,T mediated 
by the neutral weak bosons depicted Figure [RD bR. The apex NG will state 



^The same process is also mediated by the Higgs boson with an amphtude oc 
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for neutral currents. The Feynman amplitude for this process is: 

\Mr = ^ Hk'h'il - l^Hk)] Wh.igv - 9a1,Hp)] (B.9) 

where 

gv = 2sm'9w-l 9a = -\ (B.IO) 

In the ultrarelativistic target approximation we get the total cross section 

a'''' = ^{9'v + 9v9A + 9l) (B.ll) 

If the colliding particle is an anti-neutrino, that is if we consider the process 
Via ^ Vi a, the amplitude will be the same than (]B.9I) with the substitution 
qa ^ —9a which leads to 

crr = ^{9'v-9v9A + 9\) (B.12) 

B.2 Pair annihilation and creation: e+e~ ^ f f 




Figure B.2: Pair annihilation and creation. 



This process is depicted in Figure IB.2I The spin average factor X for it can be 
found in textbooks and has the form: 



X 



2mlm'j[{pi +P2) 



212 



{{PiP'i){P2p'2) + {PiP2)iP2Pi) + ml{p[p2) + m){pip2) + 2mlm}} (B.13) 
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Making the substitutions: 



PiP'i = P2P2 = E"^ — pp' cos 6 
PiP'2 = P^p'i = E"^ + pp' cos 6 



P1P2 



E'+p' 



II 7712 I /2 

P1P2 = E +p 



(P1+P2)' = ^E^ 



we finally get 



'da" 



L = m^+p = rrij + p 



12 



CM 



e^a^ \p'\ 
16^6^ 



'e^ + pV' cos^ 9 + E\m) + ml) 



(B.14) 



O'CM 



e a At: \p'\ 
16^6 ^ 



(B.15) 



In the following we analyze some special but interesting limits which have been 
used in this thesis. 



Equal masses 

In this limit we have me = mf = m^p = p' and the differential cross section 
becomes 



'da' 

dn 



2^,2 



CM 



e a 



16£;6 



[E^ + p* cos' e + 2EV 



(B.16) 



Integrating over the solid angle we get 

-{2E^ + m 



<^CM 



4 3^6 



(B.17) 



where a = |^ and s = {pi + P2Y has for this special case the value s = AE'^. In 
the ultrarelativistic limit m ~ eq. flB.17p gives 



rel 47r 2 ol 



CM 



—e a 
3 s 



(B.18) 



which is the usual cross section for lepton production at high energies - see for 
instance the process e^e~ — > 
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Necir the energy threshold 

Let us define now the variables 

mf 

y - 



E E 
which can be used to recast the total cross section as: 



O'CM 



(B.19) 



(B.20) 



Since we are assuming m/ > me, the threshold is found when E ^ rrif; we can 
have the two cases 



lini acM 



4 





TTLf = me = m 



(B.21^ 



B.3 Scattering of ordinary and mirror matter: 

e/ ^ ef 




Figure B.3: Mirror and ordinary matter scattering. 



This process is shown in Figure [K3l The spin average factor X has the form: 

X -- 



c\^ 



2mlm'i[{pi — p'l 



212 



■ {(PiP2)(p>2) + (PiP2)(P2Pi) - ml{p2p'2) - m){pip\) + 2m^m5} (B.22) 

Let us assume now me = m/ = m, which implies |p| = |p'|. Making the 
substitutions: 



PiPi 



= E — p cos 9 
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P1P2 



= + p"^ cos 9 



P1P2 
2 



bi - P2 



I(0,p-P)r = -Vsin^ 



we finally get the differential cross section: 
da 



dn 32 sin^ (f ) E V 



+ / (I + cos^ 6) + AE'^p'^ cos^ ^ ^ 1 ~ 2m'^{E'^ - cos0) + 2m^ 



;B.23) 



(B.24) 
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